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Abstract. This paper is concerned with a complete asymptotic analysis as 
!/ — > of the Munk equation dx'tp — uA'^ip = r in a domain fl C R^, supple- 
mented with homogeneous boundary conditions for ip and dnip- This equation 
is a simple model for the circulation of currents in closed basins. A crude 
analysis shows that as u —> 0, the weak limit of ip satisfies the so-called Sver- 
drup transport equation inside the domain, namely dxip'^ = t, while boundary 
layers appear in the vicinity of the boundary. 

These boundary layers, which are the main center of interest of the present 
paper, exhibit several types of peculiar behaviour. First, the size of the bound- 
ary layer on the western and eastern boundary, which had already been com- 
puted by several authors, becomes formally very large as one approaches north- 
ern and southern portions of the boudary, i.e. pieces of the boundary on which 
the normal is vertical. This phenomenon is known as geostrophic degeneracy. 
In order to avoid such singular behaviour, previous studies imposed restric- 
tive assumptions on the domain f2 and on the forcing term r. Here, we prove 
that a superposition of two boundary layers occurs in the vicinity of such 
points: the classical western or eastern boundary layers, and some northern 
or southern boundary layers, whose derivation is completely new. The size 
of northern/southern boundary layers is much larger than the one of western 
boundary layers {u^''^ vs. u^'^). We explain in detail how the superposition 
takes place, depending on the geometry of the boundary. 

Moreover, when the domain Q is not connex in the x direction, ip^ is not 
continuous in Q, and singular layers appear in order to correct its discontinu- 
ities. These singular layer are concentrated in the vicinity of horizontal lines, 
and therefore penetrate the interior of the domain Q. Hence we exhibit some 
kind of boundary layer separation. However, we emphasize that we remain 
able to prove a convergence theorem, so that the boundary layers somehow 
remain stable, in spite of the separation. 

Eventually, the effect of boundary layers is non-local in several aspects. 
On the first hand, for algebraic reasons, the boundary layer equation is radi- 
cally different on the west and east parts of the boundary. As a consequence, 
the Sverdrup equation is endowed with a Dirichlet condition on the East 
boundary, and no condition on the West boundary. Therefore western and 
eastern boundary layers have in fact an influence on the whole domain Q, and 
not only near the boundary. On the second hand, the northern and southern 
boundary layer profiles obey a propagation equation, where the space variable 
X plays the role of time, and are therefore not local. 
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CHAPTER 1 



Introduction 



The present paper is mainly concerned with mathematical methods investi- 
gating singular behaviours on the boundary of a domain Q, when the size of the 
boundary layer depends strongly on its localization, and more precisely when it 
becomes degenerate on some part of the boundary. 

Such a situation has been depicted in many former works, but never really dealt 
with insofar as additional assumptions were often made to guarantee that boundary 
terms vanish in the vicinity of the singularity. This is the case for instance of our 
study '5 of the /3-plane model for rotating fluids in a thin layer where we suppose 
that the wind forcing vanishes at the equator. The same type of assumption was 
used in the paper J18j of F. Rousset which investigates the behaviour of Ekman- 
Hartmann boundary layers on the sphere. This holds also true for the work of 
Desjardins and Grenier [6^ on Munk and Stommel layers where it is assumed that 
the Ekman pumping (which is directly related to the wind forcing) is zero in the 
vicinity of the Northern and Southern coasts. The difficulty was pointed out by D. 
Gerard- Varet and T. Paul in 8 . 

We intend here to get rid of this non physical assumption on the forcing, and to 
obtain a mathematical description of the singular boundary layers. More generally, 
we would like to understand how to capture the effects of the geometry in such 
problems of singular perturbations on domains with boundaries. 

1.1. Munk boundary layers 

The equation we will consider, the so-called Munk equation, can be written as 
follows when the domain Q is simply connected; 

Bd^^ - vA^i) = Bt on fi, 
(1.1) 

^"100 = 0, (71- VV')|an = 0. 

This model comes from large-scale oceanography and is expected to provide a good 
approximation of the oceanic vorticity potential assuming that 

• the motion of the fluid is purely two-dimensional Q, C R^ (shallow-water 
approximation) , 

• the wind forcing is integrated as a source term r g W'^'°°{il) (Ekman 
pumping), 

• the Coriolis paremeter depends linearly on latitude (betaplane approxi- 
mation) , 

• nonlinear effects are negligible, 

• the Ekman pumping at the bottom is negligible. 

When O is not simply connected, the boundary conditions are slightly different. 
We refer to paragraph 1 1. 2. 3| for more details. 
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1.1.1. Derivation of the Munk equation. 

Of course it is not completely clear whether our simplified model is really rele- 
vant from the physical point of view. Let us therefore explain briefly the derivation 
of this model (following Desjardins and Grenier j6j), and track the different simpli- 
fications that should be compared with experimental data. For further discussions 
regarding the physical relevance of such models, we refer to |15j . |16j or [9]. 

• As usual in large-scale oceanography, we start with the 3D incompressible 
Navier-Stokcs equations in a rotating frame: the velocity field u is assumed to be 
divergence- free 

and to satisfy the dynamical equation 

dtu + (u ■ V)u + 2ijj /\u = \- g -\ — J^u 

P P 

expressing the fact that the fluid evolves under the combined effects of the Coriolis 
force, the pressure, the gravity and some turbulent dissipation mechanism. The 
precise formulation of this last contribution (involving a turbulent viscosity) 

-J^u = Vh^hU + vsdssu 
P 

even commonly used by physicists, has no real justification, which is a first limita- 
tion of our study. 

These equations are set in a bounded domain 

V^{xenx n/hBixh) < xs < 0} 

and supplemented by boundary conditions. On the bottom (which is described by 
the topography Hb) and on the lateral boundaries, we assume that the fluid-solid 
interaction can be catched through a no slip condition (Dirichlet condition) : 

For the sake of simplicity, the free surface is replaced by a prescribed spherical 
boundary corresponding to the depth x^ — (rigid lid approximation), and the 
effect of the wind is modeled by some non homogeneous Navier condition 

This drastic but standard simplification is investigated for instance in |14j , but has 
no rigorous justification, which is the second weakness of the model. 

• Far from the poles and the equator, i.e. around a latitude -do G (0,7r/2), 
these equations can be rewritten in scaled cartesian-like coordinates. This involves 
many non dimensional parameters characterizing the physical properties of the flow, 
especially 

• the Rossby number which measures the relative size of the Coriolis force 

U 

2wsin^o^ 

where U and L are the typical velocity of the fluid and (horizontal) length 
of observation; 
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• the aspect ratio which characterizes the geometry of the domain 

h 

h being the typical height of the ocean, and which will be assumed to be 
of the same order as the curvature parameter L/Rq {Rq being the radius 
of the Earth) . 

The Munk equation is obtained in the fast rotating limit with thin layer approxi- 
mation, that is when e — >■ 0, with the following choice of scaling 

p = 0(e), 1^3 = 0(6) 

and Hb = —h{l + erju) ■ 

• A formal asymptotic analysis (based for instance on asymptotic expansions) 
shows that the limiting flow is purely two-dimensional: we indeed obtain at leading 
order 

(1.2) W3 = 0, 93Mh = 0, V,, •M/,-0. 

Since these conditions are not compatible with the bottom and surface boundary 
conditions, one has to introduce boundary layer corrections, referred to as Ekman 
layers, which contribute to the global energy balance via Ekman pumping. 

The dynamical equation is then obtained at next order (more or less as a 
solvability condition): since Uh is divergence- free, it is indeed completely determined 
by its vorticity: 

(1.3) V-^ • {dtu + u ■ VhU - Pyu^) = V-^ • {vhAhU - (3T + rjsu^ - au) 

where —au is the Ekman pumping associated to the energy dissipation by friction 
on the bottom, tjbu^ accounts for the effects of the topography, and ~/3T is the 
source term resulting from the wind forcing. When the domain U, is not simply 



paragraph 1.2.31 



connected, (1.3) must be supplemented with further compatibility conditions (see 



Note that this formal derivation can be justified by classical energy methods [6], 
at least for well-prepared initial data (i.e. for initial data satisfying the constraint 



equations ( 1.2 )). Starting from the limiting system ( 1.2 1( 1.3 1, one can indeed build 
a smooth approximate solution, then control its L^ distance to the solution u^ of 
the Navier-Stokes equations using some strong-weak stability principle. Note that 
one even obtains a rate of convergence. 

• The apparition of Munk boundary layers and the intensification of oceanic 
currents on western coasts we would like to describe with such a model arc physical 
phenomena which are typically linear and which result from the fact that the no-slip 
Dirichlet condition on the coasts of the basin becomes a non admissible boundary 



condition if the viscous dissipation is not a leading order term in the equation ( 1.3 1. 
More precisely, we expect to exhibit such a behaviour when the remaining 
rotating term /3V"'- • {yu^) = Puy (due to the inhomogeneities of the local rotation 
vector) is large compared to vA^-^ ■ u. In particular, we expect that boundary 
layers 

• should not depend on the topographical effects and Ekman pumping (if 
/3 is large enough) 
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• decouple from the convection (at least if boundary layers are stable in the 
sense that their kinetic energy remain small) 

• are quasi-stationary, meaning that the equation governing boundary layers 
does not involve time. 

We emphasize that our study remains entirely valid if the Coriolis factor sin-i? is 
not linearized around a given latitude 'do- In this case, in equation ( |1.1[ ), (3 is not a 
parameter, but a smooth function which remains bounded away from zero (namely 
(3o cos -(9) . 



1.1.2. Asymptotic behaviour for the Munk equation. In the time- 
dependent case, Desjardins and Grenicr have proved the following 

Theorem 1 (Desjardins-Grenier). Let fl be a smooth domain, defined by 

(1.4) n := {{X,y) e R^, Xwiv) <X < Xsiy), Vmin <y < Vmax} , 

where xe,xw G C^iVmimymax)- Assume that the wind forcing t ~ V^ • T G 
L°°([0, T], i/^(f2)) vanishes identically in the vicinity of the North and of the South 

(AO) 3A > 0, T{x,y) = Oify<y,nin + >^ory>ymax->^- 

and that it is small enough, i.e. 



(Al) 



t{x' , y)dx' 
e{v) 



<Cxv.. 



where v is the (fixed) horizontal viscosity and C\ is a universal constant depending 
only on i7 and A. 

Denote by up — V^^^ any solution to the vorticity formulation of the 2D 
Navier- Stokes- Coriolis system 

dtAi} + u ■ VAV- + Pdiij) - z^AV = /3r, 

V^ian^O, (n • VV')|ao = 0. 

Then, for all N ^ N, if s is sufficiently large, there exists 

N 
1=0 

such that 

hp - <piu~((o,T),L^(a)) < cr'^^-'^^'- 

The main order term Mq"* — W^ip^ is the weak limit of up in L^ as /3 —)> oo and 
satisfies the Sverdrup relation 

(1.5) 9,^" = T, Vf'r, = 

where T e is the East boundary 

Tb = {{x,y) G dn,x^XE{y)}- 

The boundary layers uf^ are located in a band of width {v/ji)^'^ in the vicinity of 

dn. 

The approximate solution Uapp is computed starting from asymptotic expan- 
sions in terms of the Coriolis parameter /3. 
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• At order j5 in the interior of the domain, we get the Sverdrup relation 

We can then prescribe only one boundary condition, either on the Eastern 
coast or on the Western coast. 

• Since u™* does not vanish on the boundary, we then introduce boundary 
layer corrections, which are given by the balance between pd^i''^'" and 
i/A^ijj^^. Because the space of admissible (localized) boundary correc- 
tions is of dimension 1 on the Western boundary, and of dimension 2 on 
the Eastern boundary, we can recover all the boundary conditions except 
i/iPp ~ 0, which is chosen as the boundary condition for the Sverdrup 
equation. 

The proof relies then on a standard energy method. The assumption on t guar- 
antees that V''^ is smooth, and that the coefficients arising in the definition of the 
boundary layer terms are uniformly bounded, and actually small enough to be ab- 
sorbed in the viscosity term. Note that the smallness condition (Al) is required 
only to deal with the nonlinear convection term. 

From now on, we keep only one small parameter, namely v/ /3, which we rename 



as ly. In other words, we work w^ith /3 — 1; for the linear stationary equation ( 1.1 1, 
the two visions are striclty equivalent. For the non-stationary equation, the case 
i^ = 0(1), /3 ^ 1 can be dealt with provided that we get convergence for times 
of order /3^^. The previous theorem implies in particular the following (weaker) 
statement for the linear Munk equation: 



Corollary 1. Let Q be a smooth domain defined by {1. 4). Consider some wind 
forcing r = V^ • T G L°°{[0, T], iJ^(i7)) vanishing identically in the vicinity of the 
North and of the South (i.e. satisfying assumption (AO)). Denote by Ui, — V^V'i/ 
any solution to the vorticity formulation of the 2D Stokes- Coriolis system 

dtA^ + drip - vA^ip = r, 
V'|9n=0, (n • V-0)|ao = 0. 

Then, for all N ^ N, if s is sufficiently large, there exists 

N 
i=0 

such that 

\K - <plU-((0,T).L^(O)) < C;.(^-l)/«. 

The main order term u™* — V^V-''^ ^s the weak limit of Uy in 1? as i^ -^ and 
satisfies the Sverdrup relation 

(1.6) a,^° = T, vH^^ = 

and the boundary layers u^^ have typical size v^''^ . 

Our main goal in this paper is to get rid of the assumption (AO) on the forcing 
T, and to treat more general domains fi (e.g. domains with islands, or which are 
not convex in the x direction). 
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Regarding the assumption {AO), the difficuhy is that the incompatibihty be- 
tween the formal hmit equation and the boundary conditions 

is not of the same nature depending on the part of the boundary to be considered: 
at the North and at the South, the transport term is tangential to the boundary 
and therefore is not expected to be singular even in boundary layers. 

Introducing boundary layer terms to restore the boundary conditions, we will 
find typically 

dxtpE,w - ^dlil)E,w = , 
on the lateral boundaries, and 

dxipN,s - t^dyll^N.S = , 

on the horizontal boundaries. This implies in particular that western and eastern 
boundary layers should be of size i'^''^ while northern and southern boundary layers 
should be much larger, of size v^'*. Of course there are superposition zones, which 
have to be described rather precisely if we want to get an accurate approximation. 
Of course this is the first step towards the understanding of more complex 
geometries. Additional difficulties will be discussed the next section. 

1.1.3. Stability of the stationary Munk equation. Another important 
difference with [6] comes from the fact that ( |1.1[ ) is a stationary equation, so that 
classical energy methods ( |1H 120) 1 are irrelevant. 

Denote by ^app the solution to the following approximate equation 

dxipapp - vA'^i'app ^ T + St on fl, 

i'appian = 0, {n ■ V'ijjapp)\dn = 0. 

In particular, denoting 5ip = tp — ipapp, we have 

dxSib ~ lyA'^Stp = 6t on fl, 
(1.8) 

A possible strategy to get some stability estimate is to introduce weighted 
spaces as suggested by Bresch and Colin 1 . From the identities 

o ^ Jn 

(A^<5V')e^(5V' = / (A(5V)^e^ + 2 /" e^A<5V' d^StP + [ e"" Mijj 5ip, 
m Jn Jn Jn 

together with Cauchy-Schwarz inequality, we indeed deduce that 

^ / S^^e^ + ^ / (AS^re^ < - / SrSi^e^ - 2. / ASi^dJ^e^ . 
2 Jn 2 Jq Jq Jq 

Assume that St can be decomposed into 

St — Sti + St2 , 
with Sti E L'^{il), St2 G H^^iVt). Then, using the Poincare inequality, 

^rJVe" < C{\\St^\\lA\H\\l- + \\5t.2\\h-A\^H\\l-). 
n 
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while 

< C\\Mi^\\LA\dx5nL- < C\\AS^P\\%^\\Sij\\]^,\ 
It comes finally, assuming that i^ <C 1, 
(1.9) \\5^\\l. + ^WAS^Wh < C (\\6r^\\l. + Ml^Mlzl 



n 



We then define some relevant notion of approximate solution. 

Definition 1.1.1. A function ipapp G H^JJl) is an approximate solution to 



(1.1) if it satisfies the approximate equation {1.7) for some 5t £ H (fJ) such that 



St ~ Sti + St2 with 



5ri||22 +1^ ^\\dT2\\H-2 



(1-10) lim , "^„ „ " ^ "" =0. 



Plugging (1.10 1 in (1.9) we obtain 



11^ - '^appWh + i^||A(V' - i^app)\\l2 = o(||V'app||i2 + uWA^pappWh) 

meaning that ipapp ~ ip. 

By such a method, we exhibit the dominating phenomena in terms of their 
contribution to the energy balance. In particular, for the oceanic motion, we expect 
to justify the crucial role of boundary currents as they account for a macroscopic 
part of the energy. 

1.2. Geometrical preliminaries 

The goal of this section is to state the precise assumptions regarding the domain, 
under which we are able to derive a nice approximation for the Munk problem. 

As we will see in the course of the proof, error estimates depend strongly on 
the geometry of the domain, and especially on the flatness of the boundary near 
horizontal parts. 

We will therefore consider both 

• the generic case, when the slope of the tangent vector vanishes polynomi- 
ally; 

• and a special case of flat boundary, for which the exponential decay is 
prescribed. 

1.2.1. Regularity and flatness assumptions. 

First of all we assume, without loss of generality, that fl is connecteqj We 
further require that 

the boundary is a compact C* manifold, described by a finite number of charts. 

In particular, the perimeter L is finite, and the boundary of O can be parametrized 
by the arc-length s. If dfl is connected, we will often write dQ = {{x{s), y{s)) / < 
s < L}; if dn has several connected components, such a representation holds on 



If f2 is not connected, since ( |1.1[ | is a local equation, we can perform our study on every 
connected component, and we therefore obtain a result on the whole domain Q 
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every connected component of dfl. We further choose the orientation of the arc- 
length in such a way that the local frame (i, n) with n the exterior normal is direct. 

Furthermore, by the local inversion theorem, there exists 6 > such that any 
point (x, y) e il at a distance z smaller than d from dQ has a unique projection on 
the boundary. In other words, on such a tube, the arc-length s and the distance to 
the boundary z form a nice system of coordinates. 

We thus introduce some truncation function xo G C^{[0, +oo), [0, 1]) such that 

Suppxo C [0, S) and xo = 1 on [0, S/2] . 

As mentioned in the previous heuristic study, we expect the boundary layers 
to be singular on the horizontal parts of the boundary. In order to understand the 
connection between both types of boundary layers and to get a nice approximation 
of the solution to the Munk problem, we therefore need precise information on the 
profile of the boundary near horizontal parts. 

We assume that the horizontal part F^v U F5 of the boundary consists in a finite 
number of intervals (possibly reduced to points where the tangent to the boundary 
is horizontal). Denote by d{s) the oriented angle between the horizontal vector ex 
and the exterior normal n. By definition, we set 

Fjv := {s e dn, cosd{s) = and sin 6l(s) = 1}, 

F5 := {s e dn, cos 9{s) = and sin6'(s) = -1}. 

We also introduce the following notation (see figure [lT]): let Si < S2 < • • • < Sfc 
such that 

cos(6l(s,)) = l<i<k 
and for all i G {1, • • • ,fc}, either cos{9{s)) is identically zero or cos(6'(s)) does 
not vanish between Si and s^+i. Throughout the article, we use the conventions 
Sfc+i — Si, Sq — Sfc. We further define some partition of unity {pi)i<i<k 

k 

Supppi c]sj_i,Si+i[. 
We denote by F^;, F^y the East and West boundaries of the domain: 
Fb := {.s e 917, cos6i(s)>0}, 
Tw ■■= {s (^ dVL, cos6i(s)<0}. 
Eventually, let 
(1.12) I+:^{ie {!,■■■ ,fc}, s,£dTE}. 

The profile assumption states then as follows: for any Si, for cr = ± such that 
cos9{s) ^ on [si,Siai], 
(i) either there exists n > 1 and C 7^ (depending on i and a) such that as 

C 

cos 6'(s) |(s-Si)") 

{H2{i)) ^ 

and e^^\s) - —is - sA"-^ ior 1 < I < inf(3,n); 

[n — ly. 
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Figure 1.1. Positions of the abscissa s, 



(ii) or there exists C ^ Q and a > (also depending on i and a) such that as 



S r Si^ S Kz ySi^ ^ialjj 



{H2{ti)) 



cos6l(s)-Ce-"/l''~'''l, 

a/\s — Si I 



and 9'{s) ^ —e 

[s- SiY 



The first situation corresponds to the generic case when the cancehation is of finite 
order. The second one is an example of infinite order cancellation: in that case, 
which is important since it is the archetype of C°° boundary with flat parts, we 
prescribe the exponential decay because there is no general formula for error esti- 
mates. Notice that we do not require the behaviour of 9 to be the same on both 
sides of Si, provided the function 9 belongs ioC^{d^), so that {x{s),y{s)) G C^(90). 
This profile assumption will essentially guarantee that, up to a small truncation, 
we will be able to lift boundary conditions either by East/West boundary layers, 
or by North/South boundary layers at any point of the boundary. This is therefore 
the main point to get rid of assumption ( AO) . 



1.2.2. Singularity lines. 

In the case when the domain 57 is not convex in the x direction, we will see 
that the asymptotic picture is much more complex, especially because the solution 



to the Sverdrup equation (1.5) may be strongly singular. 



We will therefore need to build singular correctors which are not localized in 
the vicinity of the boundary. This construction is rather technical, and for the 
sake of simplicity, we will use the following additional assumptions (which are quite 
general) . 
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y=yo 



The set {s, y(s)=y^ and cos e(s)=0}={t^}U{y 
is not connected. 



Figure 1.2. Example of a domain when assumption (H3) is not satisfied 



(H3) Vyo e R, the set 

{sedn, yis)=yo, cos0(s) = O} 



is a connected set (see figure |1.2[ ) . 

(H4) Let Sj G {si, • • • , s^} be a boundary point such that 

• Sj e dTEndTw; 

• and fl is not convex in a; in a neighbourhood of {x{sj),y{sj)). 
Then cos 6l(s) = 0{\s - Sj\*) for s -^ Sj. 



Assumption (H4) will be discussed in Remark 3.4.6 

We also introduce the lines S^j, across which the main order term will be 
discontinuous, which will give rise to boundary layer singularities: we set 

M 

f B := U r. 



-J' 



j=i 



where Fi, • • • , Tm are the closed connected components of Te- For j e {1, • • • , M}, 
we set 

Aj := {{x,y) e n, 3y' e R, {x',y) e T^ and {tx + (1 - t)x',y) e f7 Vi e]0, 1[} . 

We have clearly 

M 

n^lJA,. 



We also define (see Figure 1.3) 
(1.13) 



Sjj := A, n Aj for i =^ j, S := |J E, 
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Figure 1.3. Discontinuity lines of a general domain il 



It can be easily checked that every set T,ij is either empty or a horizontal line with 



ordinate yij such that there exist 



e R with {xi,yij) € Ti, (x^, y^) G Tj, and 



either {xi,yij) e dVi or {xj,ytj) G i9rj. 

Eventually, we parametrize every set Tj by a graph, namely 



r. 



{(2;i;(y),|/), yLn < 2/ < 2/max} 



1.2.3. Domains with islands. 

When the domain il is not simply connected, the boundary conditions on dfl 
are slightly different. This case has been studied in particular in [2 , where the 
authors investigate the weak limit of the Munk equation in a domain with islands. 
Let r^i, • • • flK be C^ simply connected domains of R^, such that 

• rii (^rii for i > 2; 

• n, n n^ ^ $ ioT 2 < i < j < K] 

• f] := fii \ Ui>2ili satisfies (H1)-(H4). 

Let Ci — dfli for i > 1. Notice that of course, the presence of islands gives rise to 
discontinuity lines S as described in the preceding paragraph. 
Then the Munk equation can be written as 



(1.14) 



dxi^ ~ v^ ■(/; = r in 17, 

dni) = on 917, 

V'ICi = 0, Vic. = c* for i > 2. 
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The constants Ci are different from zero in general: indeed, the condition u-n^gfi = 0, 
where u — V^?/; is the current velocity, becomes ip = constant on every connected 
component of dfl. However, the constants are not required to be all equal. In fact, 
the values of C2, ■ ■ • ,ck are dictated by compatibility conditions, namely 






(1.15) Vj>2, ly d.nM 

'Cj JCj 

where t = curl T. 

In the rest of this paragraph, we will first present in a few words where con- 



dition (1.15) comes from. Then, we will prove that the coefficients Ci are uniquely 



determined by the condition (1.15). Eventually, we will explain how the study of 
(1.14 )-( 1.15 ) can be deduced from the one of ( |1.1[ ) in an arbitrary domain satisfying 
(H1)-(H4). 



• The compatibility condition ( 1.15 ) is inherited from the Navier- Stokes system 
satisfied by m = 'V'^ip (see [T3| for a similar argument in the inviscid case). Indeed, 
we start from the stationary Stokes-Coriolis system in dimension 2, with /3-plane 
approximation, namely 

-(1 + (3ey)u^ - vAu + Vp = f, in f7, 

(1-16) divu = 0inf7, 

u\dn = 0. 

The idea is to take the curl of the first equation in order to get rid of the pressure 
term. However, for $ g L^(rj)^, the identity curl $ = does not imply the existence 
of g € H^{VL) satisfying $ = Vq. Indeed, $ is a gradient if and only if its circulation 
around any closed contour C in 51 vanishes: 

where t is the tangent vector to the curve C. If curl$ = 0, this condition becomes 
J„ $ • i = for all i > 2. Hence we obtain 



3q^H^{n), $ = Vq 



curl $ = 

and /p $ • t = for all i > 2. 



Therefore (1.16) is equivalent to 

f3u2 — I' A curlu = curlT in fl, 
divu = 0, 
u\an = 0, 



1/ / Au-t 

'Ci 



T-t^Q Vi > 2. 



This amounts to (1.14|-(1.15l with u — W-^ip (notice that, as above, the existence 
of ip is ensured by the divergence free condition on u and by the Dirichlet boundary 
conditions). 

• In order to check that the constants Ci are well-defined, we introduce, as in 
[2] , the functions ipir ' ' ■, 4'k such that 

dx4'i ~ v^ 4'! — 'T — curlT, 
V'lian = 0, dnipi\dn = 0, 



(1.17) 
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and for i> 2, 

(1.18) 



V-.ic^^ yje{lr-- ,K}. 



The existence and uniqueness of the functions ipi G H^(fl) follow for instance from 
the Lax-Milgram Lemma. We then have 



(1.19) ^P = lPl+Y,C^^^, 

i>2 

and condition ( |1.15[ ) reads, for all j > 2, 

^Cii/ / dn/^i^i = -^ / f^«^^i 

i>2 '^J '^J 



r-^-n. 



c. 



Therefore we have to prove that the matrix (v ^q. dn^4'i)2<i,j<K is invertible. We 
prove in fact that for all (02, • • • , Qk) & R^~^ 



(1.20) 



y^o-zajv / dnAtP, = -v I lAV'aP < 0, 



where ipa = X]j=2 ^iV'i- The invertibility of the matrix {v Jq. dnA^i)2<i,j<K follows 
easily. 

We have 

d^'4)a - vA'^il^a = 0, 
dn^a\dn = 0, 

Therefore, since Ci is a closed contour, 

1 ^ /■ 



/ 



4=2 



and 



Ja Jn j^2 "'C'i 

= ^ j |Af/'a|^ + y]aiaj / dnAil). 



The identity (1.20) follows. 



• The formula (1.19) shows that it is sufficient to understand the asymptotic 
behaviour of the functions ipi. Indeed, the coefficients ct (which depend on v) are 
obtained as the solutions of a linear system involving the functions ipi. It is proved 
in [2] that for any domain V C fl such that F n E = 0, where S is defined by ( 1.13 ) 

iPi -^ V° in i^(f^), 
V'l -^ V° in L^iV), 
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Figure 1.4. Definition of the domains Bi 



where ip^ is the solution of the Sverdrup equation (1.51, and 

A -^ 1b, in L^{V), 
where Bi is defined by 

B^ := {(x, y) e n, 3y' e R, (x', y) G Q n f e 

and (te + (1 - t)a;', y) enVte (0, 1)} . 

This result is in fact sufficient to compute the asymptotic limit of i^ J„ dnAtpj, 
and thereby, of the coefficients c^, which converge towards some constants q. We 
refer to section |2.5| for more details. We will go one step further in the present 
paper, since we are able to compute an asymptotic development for the functions 
ipi, and therefore give a rate of convergence for the coefficients Ci and the function 



1.3. Main approximation result 

We will prove the existence of approximate solutions in the form 



(1.21) 
where 



V'a 



V-"* + i^E,W + ^JV,S + i' 



• -0™* is a regularization of f/'" {+J2i=2'^i'^Bi '^hen the domain fl has is- 



lands), and ijj^ is as previously the solution to the Sverdrup equation ( 1.5 ) 
vanishing on the eastern boundary 

V'ite = where Te = {x e dfl / e^ ■ n > 0} 

Furthermore, Ht/^™* - <||i2 = o(||V'™*|!l2); 

• il'E,w groups together eastern and western boundary terms, the decay of 
which is of order z/^/'^; 

• ^N,s is the contribution of southern and northern boundary layers, the 
size of which is typically z/^/**; 
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• and ip is a boundary layer term, located in horizontal bands of width 
1/^'^ in the vicinity of every singular surface 'Sij. 

Notice that ^™*, ipE,w, '^n,s and ip^ do not have the same sizes in L^ and 
H^: typically, if the boundary condition to be lifted is of order 1, the size in i^ of 
a boundary layer type term is A~^/^, where A~^ is the size of the boundary layer, 
while its size in H^ is A^^^. As a consequence, we roughly expect that 

(1.22) U'"'\\l^ » Un,s\\l- ~ WWl- » Ue.wWl-, 

•^ ^ " V ' " V ' 

~1 ^1/1/8 ^1/1/6 

while 

(1.23) WWh^ ^ UnMIh- - W'Wh'^ « Ue,w\\h-, 

'^ V ' "^ V ■' 

Our result is the following; 

Theorem 2. Assume that the domain fl satisfies assum,ptions (H1)-(H4-)- Con- 
sider a non trivial forcing t = curlT G W^'^ (il) , and let ip be the solution of the 



Munk equation (1.1 1 (or ( 1.14|-( 1.15 I when fl has islands). 



Then there exists a function ipapp of the form 

i'app = V-"* + ^E,W + i^N,S + V-'^ 

such that 

U - i'appWL- = 0{V^/^), U - 4>app\\m = 0(1^-3/8)^ 

More specifically, ^jj^-pp is a good approximation of ■0, in the following sense: let 
V d ^ be a non-empty open set. Then, for a generic forcing term t, the following 
properties hold: 

• Approximation of the interior term: if V <£ fl is such that V H T, = $, 
then 

U-^°\\LHv)=o{U°\\mv)); 

• Approximation of the S boundary layers: i/y <e 51 is such that Vn'Eij 7^ 
for some i j^ j € {1, • • • , M}, then 

• Approximation of the North and South boundary layers: if V D Tn,s 7^ 
and V n Te^w = 0, t^ n S = 0, then 

\\tp - (V'" + V'JV,s)||L2(y) = o{\\iIj" + i/jn.sWl^v)), 
IIV" - '>PnM\h^{v) = o{\\iPnM\h^{v))- 

• Approximation of the West boundary layer: ifVnT\Y ^ andVnTpfg — 

0, yns = 0, ynr£ = 0, 

IIV" - i^wWH^v) = o{\\^|;w\\H^v))■ 

In the above Theorem, the term "generic forcing" is necessarily unprecise at this 
stage. It merely ensures that the terms constructed in the approximate solution are 
not identi cally z ero. We will give a more precise assumption in the next chapter (see 
Definition 2.4.1). The interior of Fjy.s, namely Tf^g, is to be understood through 



the induced topology on dfl. 



16 1. INTRODUCTION 

Remark 1.3.1. • One of the main features of our construction lies in the 
precise description of the connection between boundary layers. In particular, we 
prove that the sizes and profiles of the North/South and East/West boundary layers 
are unrelated; the transition between both types of boundary layers occurs through 
their amplitude only. Notice that this result is rather unexpected: indeed, most works 
on boundary layers (see for instance 8\) assume that an asymptotic expansion of 
the form 

V'~V'"* + ^''^(s,A(s)z) 
holds, where the size of the boundary layer, X^^, is defined on the whole boundary 
and is continuous. Here, we exhibit a different type of asymptotic expansion, which 
shows that a superposition of two types of boundary layers occurs on the transition 
zone. Also, on this transition zone, the ratio between the sizes of the two boundary 
layers is very large, so that the asymptotic expansion above cannot hold. 

• The analysis we present in this paper can be extended without difficulty to very 
general anisotropic degenerate elliptic equations, in particular to the convection- 
diffusion equation 

dx^ — vl^tp = T 

which has been studied for special domains by Eckhaus and de Jager in [7j , followed 
by Grasman in [lOj , and more recently by Jung and Temam |12j . Note that in the 
case of the convection- diffusion equation, the maximum principle (which does not 
hold anymore in our case) can be used to prove convergence in L°° . In the paper [7], 
the authors exhibit parabolic boundary layers on the North and South boundaries, 
but only treat the case when the domain Q is a rectangle. 

Note that the tools we develop in the present paper allow to consider more 
general geometries, for which 

• there is a continuous transition betiueen lateral and horizontal boundaries, 

• there are singular interfaces 'S-ij , and even islands. 

For the sake of simplicity in the presentation, we have not dealt with East 
and West corner singularities (typically the case when Q = [—x,x] x (— y, y] is a 
rectangle), but the corresponding connection between lateral and horizontal boundary 
layers is actually simpler, and can be handled exactly as the connection between the 
S and West boundary layers, which is treated in paragraph [3^ 



In view of the energy estimate (1.91, the idea is to prove that ipapp satisfies 



equation (1.7 1 with an error term St = Sti + St2 such that 
(1.24) \\ST,h^n) = o(:.i/8), \\ST2\\H-Hn) = o{i^'''^) 



Notice that the energy estimate (1.9) allows us to capture both types of boundary 
layers and the interior term: the boundary layers of size v^''^ can be seen both in 
L^ and H^ norms, while the interior term is identified through the L'^ norm only, 
and the West boundary layer through the H'^ norm only. 
We therefore use the following notion: 

Definition 1.3.2. In the rest of the paper, we say that St = Sti + St2 is an 



admissible error term if it satisfies (^1.24) 



Since the proof of Theorem [2] is very technical, we have chosen to separate as 
much as possible the construction of ipapp and the proof of convergence. Therefore 
the organization of the paper is the following: 



1.3. MAIN APPROXIMATION RESULT 17 

• In Chapter [2J we expose the main Unes of the construction of the boundary 
layer type terms, namely iI'e,Wi '4'n,s 3'iid V'^j without going into the 
technicalities. 

• In Chapter [Sj we give all the necessary details for the construction, and 
we estimate the sizes of the four terms defining ipapp- 

• Eventually, in Chapter El we prove the estimate on the error term (1.24), 
which entails that ^papp is an approximate solution. 

Moreover, as a very large number of notations are introduced throughout the paper, 
an index of notations is available after the Appendix. 



CHAPTER 2 

Multiscale analysis 

Searching as usual an approximate solution to the Munk equation ( |1.1[ ) 

dx'4' — J^A^-0 = T 
in the form 

with Tp'^ satisfying the Sverdrup relation 

we see that both ip'^ and "0^^ present singularities near the North and South bound- 
aries of the domain, i.e. as cosO vanishes, as well as on the interfaces S^j. 

• First, the main term of the approximate solution, '0'^, is singular near all 
"North-East" and "South-East corners". More precisely, combining the 
integral definition of ij:^ together with equivalents for the coordinates of 
boundary points given in Appendix A, we see that the y derivatives of -0° 
explode near such corners. 

• Moreover, as we explained in paragraph |1.1.2[ the size of the boundary 
layer becomes much larger as cos^ — > 0, going from t^^/"^ to v^^^. This 
also creates strong singularities in the boundary layer terms, which are 
completely independent from the singularity described above. In fact, 
there is a small zone in which both boundary layers coexist and are related 
to one another through their amplitude. 

• Finally, for complex domains, i.e. for domains where the closure of the 
East boundary F^ is not connected, the solution of the interior problem -0° 
has discontinuities across the horizontal lines E^. These discontinuities 
gives rise to a "boundary layer singularity", which is apparented to 
the North and South boundary terms, the size of which is therefore v^''^ . 

Because of these three types of singularities, the construction of the solution is 
quite technical. We therefore start with a brief description of all kinds of boundary 



type terms (see Figure 2.1) 



2.1. Local coordinates and the boundary layer equation 

As usual in linear singular perturbation problems, we build boundary layer 
correctors as solutions to the homogeneous linear equation 

(2.1) 9^V-i^AV = 0, 

19 



20 



2. MULTISCALE ANALYSIS 



North boundary layer, size v" 




Discontinuity 
boundary layer, 



West boundary layers, size v 



Figure 2.1. Overall view of the boundary layers 



localized in the vicinity of the boundary (therefore depending in a singular way of 
the distance z to the boundary) 



V- = v-^^ls, \z) with a,v^^ = A^zV-^^ > dsi>^^ 

i)^^ -> as Z ^ oo, 
and lifting boundary conditions 

Note that the parameter A ^ 1 is expected to measure the inverse size of the 
boundary layer, it can therefore depend on s. 

It is then natural to rewrite the homogeneous Munk equation in terms of the 
local coordinates (2:, s). We have by definition of Q, s and z, 



(2.2) 



so that 



/ dx dx\ 

dz ds 

dy^ dy 

\dz ds J 



— cost 

— sin t 



{l + ze')sine 
-(l + z0')cos( 



/ dz dz\ 

dx dy 

ds ds 

\dx dy) 



— cos y — sm f 

(1 + ^6*')"^ sin 6* -(l + z6'')-icos( 
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We therefore deduce that the jacobian of the change of variables {x,y) -^ {z,s) is 
equal to (1 + z6')~^. Furthermore, 



i.BL _ „„„/)\ a i.BL Smf 



9^V"'^ = -cos6'AazV' + 



1 + x-^ze' 



A 



On most part of the boundary the second term is negligible compared to the first 
one. Nevertheless we see immediately that on horizontal parts the first term is zero, 
so that we have to keep the second one. Since A 3> 1, we approximate the jacobian 
term (1 + X-^ZO')-^ by 1. 

Similar computations allow to express the bilaplacian in terms of the local 
coordinates (which involves more or less twenty terms). However in the boundary 
layers, we expect that the leading order term is the fourth derivative with respect 
to z 

AV^ = a49|7A^^ + 0(A3). 

We will thus consider only this term and check a posteriori that the contribution of 
other terms is indeed negligible. Note that, as we want to estimate the H~^ norm 
of the remainder, we will only need to express the laplacian in local coordinates. 

We will therefore define boundary layer correctors as (approximate) solutions 
to the equation 

(2.3) - A cos edz^'^^ + sin e{dsij^^ + ^Zdzi^^^) ~ i^X'^d^z^^^ = 0. 

A 

This equation remaining still complicated, we will actually consider two regimes 
depending on the precise localization on the boundary. For each one of these 



regimes, we will neglect one of the first two terms in (2.3) (i.e. part of the dx 
derivative), so that the size of A as well as the profile of ip"^ will be different. Of 
course we will need to check a posteriori that the term which has been neglected 
can be dealt with as a remainder in the energy estimate. 

2.2. East and West boundary layers 

In this section, we construct the boundary layers on the lateral sides of the 
domain. We retrieve rigorously the result announced in the introduction, namely 
the intensification of Western boundary currents and the dissymctry between the 
East and West coasts. 

2.2.1. The scaled equation. 

Along the East and West coasts of the domain, on intervals on which cosO 
remains bounded away from zero, it can be expected that the main terms in the 



boundary layer equation (2.3) are X cos 9dzf{s,Xz) and z/A^9^/(s, Az). Hence we 
take A such that 

n\\ 1/3 

cosoir 



(2.4) Xe,w = 

V ^ 

and / such that 

(2.5) dy^~sign{cos{0))dzf. 

We recall that Z is the rescaled boundary layer variable {Z = Xz), so that Z E 
(0,cx)). 
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We look for solutions of the above equation which decay as Z — >■ oo. Conse- 
quently, the dimension of the vector space of solutions of the simplified boundary 



layer equation (2.5) depends on the sign of cos(6'): 



If cos6' > (East coast), decaying solutions of equation (2.5) are of the 
form 

f{s,Z)^A{s)eM-Z) 



• If co&O < (West coast), decaying solutions of equation (2.5) are of the 
form 

f{s, Z) ^ A+{s) exp(-e'^/3z) + A_{s) exp(-e-"/3z). 

Notice that we retrieve the dissymetry between the East and West coasts: indeed, 
only one boundary condition can be lifted on the East boundary, whereas two 
boundary conditions (namely, the traces of ip° and dnip'^) can be lifted on the West 
boundary. As a consequence, ip'^ must vanish at first order on the East coast, so 
that the role of V''^^ on Te is merely to correct the trace of 9„'0". Note that, in 
order that the trace and the normal derivative of ■0" + ^'^^ are exactly on Te, 
we will actually need an additional corrector, which is built in the next chapter. 

In first approximation, the boundary layer terms on the East and West coasts, 
denoted respectively by ipE and 'ipWi are thus defined by 

ijE{s,Z)=A{s)exp{-Z), 

^w{s,Z) =Y,A±{s)exp{-e^'^/^Z), 
± 

where the coefficients A, A_|_ , A^ ensure that the trace and the normal derivative of 
ip'^{s, z) + 4>^^{s, Xe,wz) vanish at main order on the East and West coasts. This 
leads to 

(2.6) A{s) = A^i9„V|ao(5), 

'A+is)\ ^ ( I 1 y'f-^'^oni^T 

Let us emphasize that the precise value of A on East boundaries is in fact irrelevant 
in energy estimates, since equation (2.6) implies that A — 0(i/^") on zones where 
cos does not vanish. Therefore the East boundary layer itself is not captured 
by energy estimates. But its incidence on the interior term, through the fact that 
■0Pp = 0, is clearly seen in the L^ estimate. 

2.2.2. Domain of validity. 



(2.7) 



Simplifying ( |2.3[ ) into ( |2.5[ ), we have neglected the terms X'X^^Zdzf and dsf- 
As / is an exponential profile, the term corresponding to X' X^^Zdzf is smaller 
than vX'^d'^f and Xcos{9)f as long as 



, ^4, |cOSf?|4/3 



which leads to 

(2.8) i.i/3|0'| |cos0|-^/3<l. 
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As for the term dgf, if V'lan and Sni/'ign are smooth with respect to s, the 
corresponding error terms can be neglected as long as 

A|cos(6i)|>l, i.e. |cos(6')|>:/i/^ 



This last condition is less stringent than (2.8). Hence we only keep (2.81 in order 



to determine the interval of validity of the construction. 

The intervals on which the East and West boundary layers are defined follow 



from the validity condition (2.8l. More precisely, if cos 61 ^ on (sj,Sj+i) (West or 
East coasts), we set 

st := sup |. G (^.., fii^) , ^i/3|0'| I cos^r^/^ > l| , 

^r+i = inf {.s e (^^^^,..+1) ,-^/'|e'l \cos9\-y' > 1 

Notice that s^, s^-^ depend on v and are well-defined as long as 6* e C^([0, L]). 
By definition. 



yl/3|fl'| lr.nQfl|-7/3 ^ -I Ve^Cc+«- 



iy-'"\ti\\cosV\ ■'"<! Vs e (s+,s^_^J, 



and it is easily proved that cos(6'(sf )) vanish as i^ — ?► 0, so that 



lim s+ = s,, lims^_^i = s,;+i. 



Remark 2.2.1. Notice that (2.8) is not satisfied in the vicinity of s^ ,3^,-^^. 



Moreover, the derivatives of V'laa "■''T-d 9nip?gQ with respect to s may become very 
large as s approaches Si and Si+i. Chapter^ is devoted to the estimation of the 
corresponding error terms and to the proof of their admissibility. 

2.3. North and South boundary layers 

In this section, wc construct the boundary layer terms near the intervals where 
cos 9 vanishes. Notice that these intervals may in fact be reduced to single points. 

2.3.1. The scaled equation. We consider in this section an interval (s^, s^+i) 
on which cos 9 is identically zero and sin 9 = — 1 (South boundary), so that equation 



(2.3) becomes 

(2.9) dsf + jZdzf + i^x''dy = o. 



As in section [2^21 we have to choose A and / so that the above equation is satisfied. 
The simplest choice is to take 

(2.10) Ajv,s := ^"'/'^ 

so that A' = 0, and the boundary layer equation becomes a diffusion-like equa- 
tion, with the arc-length s playing the role of the time variable. Note that such 
degenerate parabolic boundary layers have been exhibited in [7] for instance. 
This raises several questions: 



What is the direction of propagation of "time" in (2.9)? In other words, 
what is the initial data for (2.9)? 



• Is equation (2.9) well-posed? 

• What is the domain of definition (in s) of the South boundary layer term? 

• How are the South, East and West boundary layers connected? 
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We will prove in this section that equations of the type 
5s/ + 5|/ = 0, se(0,T), Z>0, 

J\s—0 Jin 

f\z=o ~ fa, dzf\z=o = /i) 

with /in € i^(0, oo), /o,/i e W^^'°°(0,r), are well-posed, and we will give some 
energy estimates on the solutions of such equations. In the present context, this 



means that on South boundary layers, equation (2.91 is a forward equation (in s), 
while on North boundary layers it becomes a backward equation. This is consistent 
with the definition of the interior term ip'^: in all cases, the boundary condition in 
s is prescribed on the East end of the interval. 

Let us also recall that the domain of validity of the West and East boundary 
layer terms does not reach the zone where cos6' = (see (2.8)). As a consequence, 



the South boundary layer term must be defined for s < a^ < Si and s > cr^^i > Si+i 
where the arc-lengths a^ , (J^_^i will be defined later on and should satisfy 



< s7 



>fT; 



This requires to slightly modify the boundary layer equation. From now on, 
we define the South boundary layer term ipg as the solution of 

Hui u sm a 

(2.11) 



dsi^s 



-1/4 cos 61 



i^S\Z=Q = 

with "^Q, *i to be defined later on (see ( |2.14| )) 



0, s> al, Z >0, 

3L 

Q\z=a 



*o, az^gfe=n = *i 



2.3.2. Study of the boundary layer equation (2.11) 



In this paragraph, we give a well-posedness result for equation 
dJ + b{s)dzf + fiis)dy = 0, s€(0,T), Z>0, 

(2.12) f\s=0 = fin, 

f\z=o = /o, dzf\z=o — /i • 
We first note that, up to lifting the boundary conditions, we are brought back 
to the study of the same parabolic equation with homogeneous boundary conditions 



and with a source term: indeed / is a solution of (2.12) if and only if 



gis, Z) := f{s, Z) - Us){Z + 1) cxp(-Z) - /i(s)Zexp(~Z) 

is a solution of 

d,g + h{s)dzg + ti{s)d%g = S{s,Z), s e {0,T), Z > 0, 

g\s=Q = 9in, 

g\z=o = 0, dzg\z=o = 0, 
where 

g^n :- f^n - /o(0) (Z + 1) exp(-Z) - /i (0)Z CXp(-Z) 

and 

5(s,Z) := -f^{s){Z+l)cM-Z)-f[{s)ZcM-Z) 

+bis)fo{s)ZcM-Z) - b{.s),h{s){l - Z)cM-Z) 
-/x(s)/o(s)(Z - 3) exp(-Z) - nis)his){Z - 4) cxp(-Z). 
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Our precise result is then the foUowing: 

Lemma 2.3.1. Let T > 0, and let fi e L°°(0,r), b e L°°{0,T), foji e 
C\[0,T]), 5„ e i2(0,(X)), and S G L^{{0,T),L^(R+)). 
Assume that there exists fiQ > such that 

m(s) > /io for a.e. s G (0,r). 

Then the equation 

ds9 + b{s)dzg + m(s)9|.9 = S{s, Z), sG (0,T), Z > 0, 
(2.13) g|s=o = .9m, 

g\z=o = 0, dzg\z=o = 0, 

has a unique solution g G C([0, T], L^(R_|-))ni^((0, T), iJ^(R+)) which satisfies the 
energy estimate 

l\\gi-s)\\l^^.^ + l\i-sWlgi--^')\\hiTc+^ds' 



< Hn\\hf^R^)+ (^J^ \\Sis')\U^R+)ds' 



2 



Proof. Let 

t2( 



V:={ve i/nO,oo), v\z=o = 0, dzv\z=o = 0}. 
For u,v £ V, s £ (0, T), define the quadratic form 



OO /"OO 

2„,,"7^ a2. 



ais, u,v)= / 6(s)azu(^) w(^) dZ + fi{s) / a^u(Z) ^^^(Z) dZ. 
Jo -/o 

Then for all u,v £ V, the map 

s i-> a(s,u, w) 
is measurable, and for almost every s £ (0, T) 

\a{s,u,v)\ < (||fo||L=c(o,T) + llAi||L-(o,T))ll'"llH2(R,^)||i;||H2(R,^) 
a{s,u,u) > Moll9|w||i2(R,^) 

> mo(I|w|Ih2(r,^) - ||azw|ii2(R^) - ||u|ii2(R,^)). 

Notice that if u e V, 

n QQ /• GO 

I|5zt^||i2 - / {dzuf - - / u5|u 
Jo Jo 

< Il"lli2 ||g|M|ii2 

-22' 
We infer eventually that 

/ \ ^ ^0 II l|2 3/io II ||2 

a(s,u,u) > y||u||ff2(R.^) - ^||w|Il2(r+)- 

Using Theorem 10.9 by J.L. Lions in J3], we infer that there exists a unique so- 
lution g e L2((0,r),y) nC([0,T],L2(R^)) of equation ^A^ such that d^g € 
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Multiplying equation ( |2.13[ ) by g and integrating on R 
1 d 



2ds 



OO /'OO /'OO 

-/o Jo 



we infer 

^(s, Z)g{s, Z) dZ. 



Therefore 



d 
ds 



Ms)\\l^- 



^'^{s')\\^l9{s')\\hds' 



< 
< 



ll^(^)IU: 

V2\\S{s)U2 



1 



9i^)\\h 



K^')\\dl9{s')\\hds' 



1/2 



Integrating with respect to s leads then to the desired inequality. 



D 



We deduce easily that (2.11 1 has a unique solution on any interval of the form 
((Tj^, cTj^]^) such that sin0 does not vanish on [Cj^, cr,^j^], and that it satisfies 

IIV'5(s)|U^(R+) < ^(1*0(^)1 + 1*1(5)1 + 1^^0(^,^)1 + 1*1^^)1) 



+ C 

+c 



i\^',is')\ + \^[{s')\)ds' 



with ^(s) == -l/sin6'(s), b{s) = -j/^i/"* cos6'(s)/ sin6'(s). 

2.3.3. Boundary conditions for s e {si, Si+i). In order to satisfy the bound- 
ary conditions 

■0|ao = 0, dntp\dn = 0, 
the South boundary layer term constructed above must be such that for s e 

(Si,S,+l), 

iy-'/^dz^s{s,Z^O)^-d^yj'^s^{s). 



In equation (2.11), we therefore take, for s G (si,Si+i), 

(2.14) ^o{s) := -^lonis): 

(2.15) ^,{s):=i^'^%^^ani^). 

There remains to define *o ^nd ^I^i for s e {<T^,Si] U [5^+1,0-^^-^). 

2.3.4. Connection with East and West boundary layers. 

• In the most simple cases, the North or South boundaries are connected on 
the one hand to the East boundary, and on the other hand to the West boundary. 
This corresponds to the situation when 

- the North boundary y — jji is a local maximum of the ordinate; 

- the South boundary y = y^ is a local minimum of the ordinate. 

Without loss of generality, we assume that the corresponding piece of di^ is a 
"South" boundary, meaning that cos6' < in a neighbourhood on the right of Si+i, 
cos 6 > in a neighbourhood on the right of Si, and sin6' = — 1 for s S [s^, s^+i]. 

The connection with the East boundary is fairly simple: we will introduce 
some truncation Xv of r close to the East corner (parametrized by Si). Therefore 
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the solution to the transport equation dxtpt — '''Xv is identically zero in a vicinity 
of {xi.yi), so that we merely require V's|s=ct" ~ *-*■ 

Concerning the connection with the West boundary layer, the situation is not 
as straightforward. In order to satisfy the boundary conditions 

V'lao = 0, dn-4^\dQ. = 0, 

the West and South boundary layer terms must be such that 

i^wis, Z^0)+ ^5(s, Z = Q) = -V?|an(s), 

\wdz^w{s. Z = 0) + Xsdz^s{s, Z = 0) = -9„7/.,%^(s) 

in a vicinity of s = s^+i. As a consequence, we take 

il)w{s, Z = 0) = -(/3,+i(s)V'°|afj(s), 

dz^w{s,Z = 0) = -A^yVi+i(s)5n'0t%n(s)> 

and 

i^s{s, Z = 0) = -(1 - ip,+,{s))^%^(s), 

dzM-^,Z^ 0) - -A^^ll - ^,+i(s))a„V?|an(s), 

where V'i+i is a truncation function which we will define precisely in the next 
chapter. We emphasize that with this definition, the South and West boundary 
layers are related via their amplitude only: in particular, the sizes and profiles of 
the boundary layers are not related. 

• In more complex cases, the line y = Hi intersects the interior of fi, and we get a 
singularity which is not only localized in the vicinity of 951. Techniques of boundary 
layers allow however to understand the qualitative behaviour of the solution in the 
vicinity of the line of singularity. This heuristic approach is presented in the next 
paragraph. 

2.4. Discontinuity zones 

Discontinuity zones occur when the East boundary Te is not connected. At 



leading order, we expect the solution to (1.1 1 to be approximated by the solution 



ip^ to the transport equation with suitable truncations Xv near East corners 

i^'tw. = 0. 

If the domain O has islands, we add to tp^ the quantity J2j=2^j^Bj- We will 
give more details on this case in the next paragraph. Note that, because of the 
truncation x^,, the main order approximation ip^ now depends (weakly) on i/. 
Therefore on every set Ai, tp^ takes the form 

'4'ti^^y) = ~ Tx„{x',y)dx'. 



Unlike the main order term in Theoremlll ip^ does not belong to H^{il) in general. It 
is indeed obvious that ■0° and dytp^ may be discontinuous across every (nonempty) 
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line Sij. More precisely, the jump of ^/"^ across a given line S^- takes the following 
form (up to an inversion of the indexes i and j) 

i'^th^j = {'rXu){x' ,y^J) dx'. 

J Xi 

The existence of such discontinuities in the main interior term is a serious im- 



pediment to energy methods. Indeed, inequality (1.9), for instance, requires the 
approximate solution to be at least in 77^(51). In fact, we will prove that this 
discontinuity gives rise to a "boundary layer singularity": the corresponding 
corrector is a boundary layer term located in the vicinity of S^ . Since the normal 
vector to 1]^- is parallel to Cy, this boundary layer term is apparented to the North 
and South boundary terms constructed in Section [2.3[ and therefore the size of the 
boundary layer is i^^'*. 

We now define what we have called in Theorem [2] a "generic forcing r" : 

Definition 2.4.1. We say that the forcing r is generic if all the following 
conditions are satisfied: 

• For all i ^ j such that S^ 7^ 0, there exists cs > such that 



for all V sufficiently small, 



^l 



K 

E 

1=2 



qIb, 



>cs; 



• // Fat, 5 =/= (for the induced topology on dfl), then for any V C Fjv.s, 
there exists cy > such that 

for all V sufficiently small, llV't'igolU^cy) ^ cy; 

• For any nonempty V C Tw, there exists cy > such that 



for all V sufficiently small. 




qIs, 



\dn 



> Cy. 



L^{V) 



• From now on, for the sake of simplicity, we restrict the presentation to the case 
when M = 2, i.e. Te has two connected components, and -0° has exactly one 
line of discontinuity. Of course, our construction can be immediately generalized 
to the case when there are more than two connected components, and therefore 
several lines of discontinuity: we merely add up the local correctors constructed in 
the vicinity of every S^j. In particular, when O has islands, there are always at 
least two discontinuity lines. We will sketch the necessary adaptations in the next 
paragraph and leave the details to the reader. 



Thus we henceforth assume that Vt has the following form (see Figure 2.2 1 



o^r^+uo^us. 



where 



• ri* are non empty, open and convex in x; 

• Vt+ = {{x,y) e 17, y > yi ov X > xi}; 

• Q.~ = {(x, y) G ri, y < yi and x < xi}; 

• E = {{x,yi) efl,x< xi} 
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Figure 2.2. Two possible configurations for the domain O 



and {xi,yi) = {x{si),y{si)) for some Si such that cos9{si) = 0. Without any 
foss of generahty, we further assume that sin0(si) = —1. This type of domain can 
correspond to two different configurations (see Figure [2?2| . 
We will use the following notations 



X 



^x<xii-y<yn X t X 



which somehow stand for lfi± : their role is to avoid any artificial singularity on 
dfl, and F^ = dfl^ n F^. We parametrize each set F^ by a graph {x^{y), y). 

As explained above, the function -0^ and its y derivative are discontinuous 
across S. More precisely, for x < xi, 

{TXiy){x',yi)dx' + / {TXu)ix',yi)dx' 

J X 

{'rXy){x' ,yi) dx' . 

'^EiVl) 

Notice in particular that the jump is constant along S. In a similar way, since 
{TXu)ix]^{y) , y) vanishes in a neighbourhood on the left of ?/i, 

[dytpt]\six,yi) = -(a;^)'(yi)(rx,.)(a;B(yi),2/i) 

rx%{yi) 

dy{TXv){x' ,yi) dx' . 
\yi) 
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We now construct the boundary layer type correctors tp , ip . The role of ip 



is to counterbalance the jump of "0^, and of its normal derivative dyip^, across S: 



When we lift these boundary conditions, we introduce a source term in the equation. 
This source term is then handled by a boundary layer type term ip^ , which has no 
discontinuity across S: ^^ S H^{U). 

At first sight, if we consider the whole singular corrector 1/)''^'+?/;^, this problem 
could seem underdetermined. Indeed there are two jump conditions to be satisfied, 
and possibly two boundary layers (on fi"*" and on fl~), each of which is the solution 



of an equation of the type (2.11), and therefore having each two entries (\E'o and 

However, when looking at the error terms, we see that the traces of dytp]j^ and 
9^V|s* appear in the energy estimates. Therefore we have further to request that 

Because of this constraint, the energy of the boundary layer discontinuity is dis- 
tributed on both sides of E. The boundary term i/)''^' + ip'^ is unequivocally de- 
fined, since roughly speaking, there are four jump conditions and four entries for 
the boundary layer terms. 

If we could have nice formulations for the North and South boundary layer 
operators (i.e. a nice formulation for the boundary layers in terms of ^1 and 
^0): we would expect to get (0^'^* -I- ip^) by a rather simple process, obtained by 
application of the previous results. Nevertheless, as 

• the boundary layer equation is highly nonlocal (propagation with respect 
to s), 

• the fictitious boundary is not perfectly horizontal, 

we are not able to proceed in such a way. We will now explain in more details our 
construction. 

2.4.1. Lifting the discontinuity. 

Let us start by introducing some notation. Consider the closed set {s G 
dfl, y{s) — yi}. We denote by /i its connected component containing si. Then /i 
is a closed interval (possibly reduced to a single point). Without loss of generality, 
we assume that si = sup/i, which means (recall that sin0(si) = —1) that si is the 
West end of /i. Hence we have either /i = {si}, or /i ~ [sq, si] with sq < si. As 



in equation (2.11 1, we introduce a point a which will be the initial point for 0''^* 
and f/;^, and which we will define more precisely in the next chapter (see (3.37)). 



In order to keep the construction as simple as possible, we use that the can- 



cellation of cos 6* near s = inf /i is strong (assumption (H4) in section 1.2), which 
implies in particular that the boundary of il^ near the junction point s — inf /i 
has C^ regularity. 

The role of V'''^' is to lift both the jumps of V'? and dy^p^ across S, and the 
traces of tp° and 9n'0? on the portion of dV, between a~ and si, so that tjj^ + 0''^* 
belongs to H'^{n). 

We have indeed the following 
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Proposition 2.4.2. Let ip be a function on fl such that 

^|n+ G H^n+), V|n- e H\n-), and [V'ls - [dn^^h = . 
Then ip belongs to H^{fl). 

Proof. The squared L^ norm being additive, the only point to be checked is 
that 

AV' = Ai/'|a+ + AV'in- 
in the sense of distributions, which is obtained by a simple duality argument. 

For any (j) £ C^(i7), because of the jump conditions on dyip and ip, we have 

JT. Jn- 

0A(^X+) + / 05j,V|s+ - / ^M^X-) - f ^dy^is- 

(/.(A(^X+) + A(Vx")) 
n 

from which we deduce that Aip e i^(il), and A?/; = A(f/;x+) + A('0x„). D 

We therefore seek ?/j''^* in the form 
(2.16) V"^=^, V) = X+X {^-^\ Hx) + b{x){y - y^)] , 



1/ 



1/4 



where the truncation x G C^ (R) is such that 

(2.17) SuppxC[-l,l], x(?/) = 1 for |y| < 1/2, 

and satisfying the conditions 

(«) V'lan+nao = -(1 - 'P)^t\dn^ for « ^ (o"", si), 
^"V'lan+nan = -(1 - 'P)dni't\dn for « ^ (cr", si), 



As in section [2. 3[ (/? is a truncation function in s, whose role is to ensure a smooth 
transition between the East /West boundary layer on [cr~,inf/i] and the disconti- 
nuity boundary layer. 

A straightforward computation then provides the following formulas for the 
coefhcients a and b 

b{x) = {l-ip{s)) -sm9{s)^n'lJJtlgais)+cos9{s)^s^Pt\^nis) 

-cos6i(s)(^Vt°|ao('S), 
a{x) = -(l-(p(s))V't°|9o(s)-5(a;(s))(y(s)-yi), 

for xi < X ^ x{s) < x(ct^), and by 

b{x) = -[dyi^^h "^ - '■ 
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For X > x{a ), we merely take a(x) = b{x) ~ 0. 



Of course, however, ip is not a solution of (1.1), even in the approximate 
sense. More precisely, we have the following lemma : 

Lemma 2.4.3. With the previous definition and notations, 

-lift ( ^ y~yi 



whe 



and 



{d, - .A'M + V'"^*) = r + Sr'^^^ ^x, ^ j + r^,, + rf,,,, 

||riJU.(n)=o(^V8)^ llriifJ|H-(o)=o(^^/«), 

St^'^\x, Y) = x^{x, yi + v^l^Y)x {y) k(a;) + v^l^}J(x)Y 

-X+(a;, yi + v^l^Y)x^'^'^ (F) \a{x) + v^l%(x)Y 
-\X^(x,yx^v^l^Y)v^l^X^^\y)h(x)- 

Definition 2.4.4. In the rest of the paper, we set 
(2.18) V™* :=V^t° + V'"^'; 

which is consistent with the statement of Theorem \^ in which ?/;*"* was described 
as an H^ regularization of ip^ . 



Sketch of proof. By Proposition 2.4.2 Vt + v belongs to H^{fl) and 



9,(^J' + V"") = Tx. + x^x 



y-vi 



{a'{x) + b\x){y~y,)), 



A(^o + ^"f*) = x+(A(V? + V'"'*))+X-(AV'?). 
Moreover, as we have the following trace identities on 9fl+ 
(A^"f*)|s = 0, (a,AV'"")|s = 0, 
applying again Proposition |2. 4. 2] we get 

^A(X+AV'*) = ^X+x(^^)(a^''(^) + 6^'H^)(2/-yi)) 

+4.^/^X+x' (^) b"{x) 

+2^'/\+x" (^) {a"{x) + b"{x){y - y,)) 

+xV^ (^) K^) + &(^)(2/-yi)] 
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We thus define 



n\,, = ^v^l\H ('^) h'\x) + 2.'/\+x" (^) («"(-) + b"{x){y - y,)) 

4,, = rix. - 1) + ux-^x ( V7?^) («*'^(^) + &*''(^)(y - V^)) 

-z.A2(x+V^0)--z.A2(x-V?)- 
so that 

• The remainder terms will be dealt with in the last chapter, using some 



technical estimates on a, b and their derivatives (see Lemmas 3.1.1 and 



(2.19) ^^ 



B.l and paragraph 4.1.2). 

• The source term (5r"" is not an admissible error term in the sense of 
Definition |1.3.2| Therefore it has to be corrected by yet another boundary 
layer type term. 

D 

2.4.2. The interior singular layer. The lifting term V'''^' has introduced 
a remainder (5r''^' which we now treat as a source term. By analogy with North 
and South boundary layers, we therefore define ip^ as the solution of the following 
equation 

d^i)'^ - dlri,'^ = -5t^'^\ x<x{a-), Y>Y_{x) 

where the function Y^ parametrizes the lower boundary of the interior singular 
layer 

• which has typical size 5y, where 5y is any parameter such that 5y 3> i^^''*, 

• which is of course included in 17. 

Note that we use here cartesian coordinates x and Y — {y — y\)lv^l^^ since the 
discontinuity line E is essentially horizontal (which is not the case for North and 
South boundary layers which are extended on macroscopic parts of the adjacent 
East and West boundaries). 

More precisely, define <^'^ and t\ by 

q+=inf{s>si, |2;(s)-y(si)| =(5J, i+ = inf{s > si, \y{s) - y{s{)\=2&y\. 

Denote by y{x) the ordinate of the point of 912 with abscissa x in the vicinity of 
{x\^y\) and of [x(^sq) ^yiySo)) . Then Y- is any C^ function satisfying the following 
conditions: 

Y^{x) = y-^^^^ forxe(x(q+),x(a-)), 

(2-20) Y_[x)=Y_>6yV-^l^ forx<x(t+), 

Vfc < 4,Vp > 1, 3C > 0, ||F-||w'=.p(s,.(.-)) < C||r-||„..,p(,(,+),,(,-)), 

where x is an arbitrary fixed point such that the vertical line x = x does not 
intersect H. 
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5Q 



x=x 



Figure 2.3. The interior singular layer 



Then the solution ip of equation (2.191 can be expressed as 



where (p^ is the solution of the backward equation 



(2.21) 



9|0^ 



= -(5T"f*(a;,y_(a;) + Z), a; < a;(cr-), 



Z>0, 



which is precisely an equation of the type (2.131. In the equation above, we still 



denote the rescaled boundary layer coordinate by Z, in order to be consistent with 



(2.13). However, in the present context, Z is just a scaled cartesian coordinate. 

Notice that ip^ and 0^ depend in fact on ly through the function y_ and the 
coefficients a and b. The asymptotic value of ip'^ as i^ vanishes will be given later 
on. 



This indeed requires a refined version of Lemma [2 . 3 . 1 1 giving high order a priori 



estimates on the solutions to the singular layer equation (2.191, as well as Lemma 



B.l controlling the sizes of a, & and their derivatives. 



2.5. The case of islands 
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When O is not simply connected, the solution of ( 1.14 )-( 1.15) is given by 



K 



where ipi is the solution of ( 1.1 1, ipi solves, ioY 2 <i < K, 

and the constants Ci satisfy a linear system whose coefficients depend on the ^^ 



Therefore, as explained in the introduction (paragraph 1.2.3) the main issue 



is to compute the asymptotic behaviour of every function -0^. Following the steps 
described in the preceding sections (and which will be developed in the next chap- 
ter), we are able to construct a function ilf^^^ for every \<i<K. Of course, there 
are a few minor changes, due to the non homogeneous boundary condition, but we 
leave those to the reader, as they do not bear any additional difficulty. 

Let us admit for the time being that Theorem [2] holds for the functions ^/^i, 
namely that we are able to construct i/'f''^ such that 

(2.22) 110, - 0^ lU. = o{v^l\ 110, - 0n|^. = o{v-^l^). 

Hence, in order to prove Theorem [2] for the function 0, it suffices to define 

0app = 0r+Ecrv'r 

i=2 



where the constants 



app 



are such that 



(2.23) |q - c'^PP\ = o(i/^/»), c^P - 0(1). 

According to (ll]22i Osl) , we have ||0r^||L2 = 0(1), 110^^11^2 = 0{iy-^/^), so 



that (2.22) and (2.231 imply 



11^ - i'apph- = o(i^^/'), 110 - ^JappWn- = oiiy-^/""). 



We now turn to the proof of (2.23), which relies on the following Lemma: 

Lemma 2.5.1. For all i G {2, • • • ,K}, let gi e C°°{Cl) such that gi = 1 in a 
neighbourhood of Ci, and Suppg^ C\Cj —% for i ^ j. 
Let us define the following matrices: 



M" := 

M^PP := 
together with the vectors 

D" := \ -V 



V I (9„A0j 

Ci 



■^]^^d^g, 



2<i,j<K 



2<i,j<K 



9„A0i 



lyo-pp ._ 



{t9o +ih^^d^gj) + 



T^ 



T' 



2<j<K 



2<j<K 
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Then the following facts hold: 

• IM"^ - M'^PP] = o(:/i/8), ID"" - D^PP] = o{v^/^); 

• M", M^PP are invertible matrices, with {M")-^ = 0(1), (M"pp)-i = 

0(i); 

• D'^PP = 0{l). 



Before addressing the proof of the lemma, let us explain why ( 2.23 1 follows: by 
definition, 

/.A 

M"" \ : = NFc = D", 

so that 

M''{c - c^'PP) = [M^PP ~ M^y^PP + 0" - D^PP, 
where c°'PP is defined by 



Using the last two items of Lemma 2.5.1 we infer that c°'Pp — 0(1). Therefore 

A/r(c-c"pp) = o(z^i/*). 



Since {M'')-^ = 0(1), ([2^ holds. 

Proof of Lemma [2.5.11 • The first step is to use the equation satisfied by 
every ipi in order to express the coefficients of iW, D'' as integrals on fl. More 
precisely, we have, since Vgi = on a neighbourhood of dfl, 



JCi Jdfl 



= 1/ / giA'^ipj + iy I VAV'j • Vg,; 



9idx^j -V I -tpjA^gi 



^jd^9i + Sij e^-n-iy ipjA g, 
n Jci Jn 



Since Ci is a closed curve, we obtain eventually 



^jdxgi -v ipjA'^g 



In a similar way. 



ly / On Alp J 
JCi 



^v I dnAipi = / {rgj + ipid^gj) + v I ViA gi 
JCi Jn Jo 



The energy estimate (1.9 1 shows that 11-0^11^2 = 0(1) for I < j < K, so that 

M" = (- 

(2.24) 



'O / 2<i.j<K 



D" = 



[rg, + i^Agj) + T^ ■ n \ + 0(zy). 

" •^'^J / 2<j<K 



Replacing every function ipj by ipj and using the error estimate (2.22), we obtain 
the first point of the lemma. 
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• To prove the second point of the Lemma, we use identity (1.201 in the intro- 
duction, from which we deduce the following coercivity inequality: there exists a 
constant C > 0, independent of i^, such that 



(2.25) 



Va = (a2, • • • , flA') e K'^^-K J2 M^ja^aj < ~C\a\' 



Of course this entails easily that (Af^)^^ is bounded. Using the first item of the 
Lemma, we infer that M°'Pp satisfies the same type of coercivity property, and 
therefore (M^^^)"^ is bounded as well. 

Let us now prove (2.25). Let a g R-^~^ be a fixed vector. Let y C 51 be an 
open set such that ^ is a neighbourhood of Fiy, namely 

• dV n dH, C Fn/ (notice that this implies that V does not intcrsct r^v.s, 
iior Te); 

• y n E = (y does not meet the discontinuity boundary layer) ; 

• Vn B, ^ d), and dV nB,nTw 7^ ^ for all 2 < i < K {V intersects the 
West boundary of every set Bi). 

Then we may write, with tp^^^ = ^ aiipil^^ , 



JQ. Jn 

Jv 



'V 

,^2 -JvnB, 



On every set V (1 Bi, by definitiorQ 



Consequently, at the main order, 



a-t ( l + Xo(z)^ 






exp 



{-Xwe^"^'z) 



IvnBi 



app 1 2 



^W 



VnBi 



Xo(^) E ^—^ exp(-AH/e±^-/3z) 



Vs 



ds dz 



> Ci^a^, / Xl,{s)ds 

JvnBinVw 

> Cal 



Gathering all the terms, wc obtain (2.25). 

• The fact that D^pp = 0(1) is obvious and follows from (1.9), which yields 

For the sake of completeness, we also give the system satisfied by c = lim^_i.o c = 
lim„^o c°-PP: 



Notice that on Bi, the main order of ^°^^^ is Is^, whose normal derivative on 9^2 is zero 
(away from the end point of S). As a consequence, several corrector terms, whose construction is 
necessary in the general case, are identically zero here. 
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2. MULTISCALE ANALYSIS 



Corollary 2. Let c be the solution of 

Mc = D, 
where 

M := i 



OBinc, 



D := 



n JCi 



2<iJ<K 



2<]<K 



Then c ~ liiny_j.o Ci • 



Proof. We pass to the limit in the formulas (2.24). It is proved in [2], and it 
can be easily deduced from our construction, that 

V-j ^ 1b, m ^2(0) for 2 < j < K, 

and therefore 

- / ijjdxgi ^ - / d^gi = - / e^; ■ n. 

Jn J Bj JdAjDCi 

Hence the identities (2.241 imply easily that 

M" -^ M, D" -^ D. 



Passing to the limit in inequality (2.25), we infer that M is also coercive, and 

D 



therefore invertible. Thus c is well-defined. The corollary follows. 



CHAPTER 3 



Construction of the approximate solution 



The approximate solutions to the Munk equation (1.1) are obtained by gath- 
ering together the different elementary pieces described in the previous chapter, 
namely 

• the interior term which is essentially the solution to the transport equation 



(1.5), regularized in the vicinity of "East corners" and of the interface S; 

• East and West boundary layer terms, which lift locally the boundary 
conditions but become singular in the vicinity of the points s,; ; 

• North and South boundary layer terms, which lift the boundary conditions 
on the horizontal parts of the boundary but in a non local way. 

We have now to understand the interplay between those different elementary pieces. 



Of course, the equation ( 1.1 ) being linear, the errors induced by all these terms 
are simply added, so that the control on the remainders in the approximate equation 
will be rather simple to obtain. The point to be stressed is that we need that each 
elementary term to be smooth enough (namely H^), with suitable controls on the 
corresponding derivatives. 



More precisely, we will consider here only one basic problem among (1.17)(1.18), 



say the case (1.17) with homogeneous boundary conditions V'|si = dntp\n = and 
forcing r. We will further assume, without loss of generality, that we have the 
simple geometry described in the previous chapter 

o = r2+uO"ui], 

where 

• H, are non empty, open and convex in x; 

• ri+ = {{x^y) G ^, y > yi or X > xi}; 

• fl^ = {(x, y) G fl, y < yi and x < xi}; 

• E = {{x,yi) eft, x < xi} 

and {xi,yi) = {x{si),y{si)) for some si such that cos9{si) — and sin0(si) = —1. 
We emphasize that the only simplification here regards notations, and that more 
complex domains satisfying assumptions {HI) — {HA) are handled exactly in the 
same way. 

In the present chapter, we will actually focus on the following difficulties 



the precise regularization process for the interior term in fl (section 3.1 ); 
the fact that the East boundary layer does not lift simultaneously both 
boundary conditions, which implies that one has to introduce an addi- 
tional corrector defined on a macroscopic domain in the vicinity of the 
East boundary but far from corners {xi,yi) where it would be singular 



(section 3.2); 
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• the connection with North/South boundary layers (even if the correspond- 
ing horizontal part of the domain is reduced to one point {xi,yi)): we in- 
deed expect these boundary layers to lift the boundary conditions both on 
the horizontal parts and on the East /West boundaries close to the corners 



{xi,yi) (section 3.3). Note that, since these North/South boundary layers 
are defined by a non local equation, we need to check that they do not 
carry any more energy beyond some point, so that they can be truncated 
and considered as local contributions; 
the truncation of the surface layer near its West end, in order that its 



trace on the West boundary is not too singular (see 3.4); 



the precise definition of the West boundary layers, which have to lift the 



traces of all previous contributions (see 3.5). 

We will also derive estimates on the terms constructed at every step. 

Note that, for boundary layers, the construction must be performed in 
a precise order, starting with East boundaries, then defining North/South 
(and surface) boundary layers, and finally lifting the West boundary con- 
ditions. This dissymmetry between East and West boundaries is similar to what 
happens at the macroscopic level for the interior term. 

3.1. The interior term 

As in the previous chapter, we define 

(3.1) ij°{x,y)^- T{x',y)dx', 



where xsiy) is the abscissa (or longitude) of some point in Te with ordinate (or 
latitude) y. We have seen that tp''^ does not belong to H^{fl) in general: indeed, 

• The function x^ has singularities near the points yi,i e /_|_; 

• Since xe takes different values on fi^ and fl^ , the function tp^ and its y 
derivative are discontinuous across E. 

Therefore, tp*^ cannot be used as such in the definition of the approximate solution 



i^^PP of (1.1) 



Let us first consider separately the domains fl^ . We remove the singularities 
by truncating the function r near the end points of T^ (with abscissa Si^i G /+). 

The size of the truncation depends on the rate of cancellation of cos 9 near Si. 
Since the rate of cancellation may be different on the left and on the right of s^, we 
seek for truncation functions with different behaviours on the left and on the right 
of every cancellation point. More precisely, let 

(3.2) x„(x, y) := Y[ Xti^ - a;j, y - y^) 

where {xi,yi) are the coordinates of the point of 90 with arc length Si. Each 
function Xi takes the following form: 
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and the function x is defined by 

1-X{x,y) :=(1-Xi(a;))(l-X2(y)) (x, y) £ R^ 
with xi,X2eC°°(R), and 

Xi(a;) = 1 if |x| > 1 + -^, xi(a;) = if \x\ < 1, 

My) = 1 if l2/l > 1, X2{y) = if \y\ < 1/2. 

The above definition calls for a few remarks: 

• Notice that the rate of cancellation in y, namely Sy, is the same for every 
function Xi- It is far from obvious a priori that such a choice can lead 
to a suitable truncation. However, in order not to further burden the 
notations, we have chosen to anticipate on this result, which follows from 
the proof of Lemma [3. 1.1| below. 

• The functions Xi all belong to C^(R^). Indeed, it can be easily checked, 
with the definition of x, that all the x derivatives of x vanish for a; = 0. 

• The X derivative of the function x is unbounded (it is of order | In 5y | in 
L°°). This choice is mandatory in the case of an exponential cancellation 
of COS0 around Si (assumption (H2ii)). If the cancellation around s^ is 
algebraic (assumption (II2i)), the function x can be any smooth function 
vanishing near zero and such that 1 — x has compact support. 

• The function x satisfies obviously 

0<X<1, 
Xix,y) = if \y\ < - and \x\ < 1, 



X(x,y) = 1 if I2/I > 1 or |a;| > 1 



1 



ln<5,| 



We set (see Figure 3.1) 



(3.3) ip°{x,y) -.^ - I Xv{x' ,y)T{x' ,y) dx' 



With this definition, ?/;? e H^{^^), but V? i H^{n): indeed 7/.° a nd dy ip^ are 
still discontinuous across S. The term ?/)''^* constructed in paragraph 3.4.1 of the 
present chapter will lift this discontinuity. 

Moreover, by definition of ip^ , we have 

with 

(3.4) Jr = r(x.-1)-^A2V°. 

Let us now give the definition of the truncation rates. We take 

(3.5) Sy -.^ v^/^\lniy\^/^\ii\\iiiy\)-^ , 

for some arbitrary exponent /3 > 0. The definition of S^ is a little more involved. 
If XE{y) is defined on both sides of yi, i.e. if Si belongs to the interior of r^, then 
we set 

(3.6) 51^^ := \xE{yi T5y) - Xi\, (5^_ := \xE{yi±5y) - x^\ where ± = s\TL0{si). 
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Figure 3.1. The function r is truncated in the hatched zones. In 
the picture, /+ = {1, 2, 3, 4}. 

If XE{y) is defined only on one side of j/i, say for y > yi, then we take 

(3.7) 6l.^_^_ = 6l._ = \xE{yi + 6y)-Xi\. 



The following result holds: 

Lemma 3.1.1. Assume that fl satisfies assumptions (HI), (H2). Let 13 > Q he 
arbitrary, and let 5y, 5x as in (3.51, (3.6 1 -(3.7 1. Then 



\\T{x.-l)\\H-Hn)=o{v^/''), 

This lemma will be proved in Chapter [4] 

We conclude this paragraph by giving some estimates on the trace of f/;^ and 
dn'4't- These estimates will be useful when we construct the boundary layer terms 
lifting these boundary conditions. For the sake of readability, we introduce the 
following majorizing function 



M{y) = l + 



1 



inf 



ie/4 



r 11 '^Sy<2\y-yi\<l 

{y-yi)ln\y-yi\\ i^i^ 



Lemma 3.1.2. Let ip'^ be defined by 



V'?(a;,y) 



2;E(y) 



{TXi^){x\y)dx'. 
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• Trace estimates on dfl: There exists a constant C , depending only on 57 and 
T, such that for all s G dO., 

\\i't\dn\\L'=-(dn) < C, 

\dni^t\dnis)\ < C I ^^n^^^i Y. My(s)-v.\<Sy + M{y{s)) 

Moreover, for all s G 957 such that {x{s),y{s)) ^ E, 
l^sV'°|an(s)l + Sy\dsdnip%n{s)\ 

- ^ VTY^ ^ hyis)-y.\<8^ + \cose{s)\M{y{s)) 
y ^' ^' iei+ 

and 

\dl4>%n{^)\+5y\dld^rt\o,M\ < ^(^ + i^^ 

"y \ "y "v 



Eventually, on the East coast, we have the following more precise estimates: 

\0'is)\' 



'^^^^"^°i-(^)' ^ ^sJ^£kwU" 



\y(s)-y,\>Sy/2- 



• Jump estimate on S; The jumps [V't'ls; [9j,^"]e o-re constant along S. More- 
over, there exists a constant C such that 

Uh\<c, |M.|<^. 

Proof. The L°° bound on V't'ign is obvious: we merely observe that 

llV't°|aollL~(cio) < IklU-ai) sup |x-x'|. 

(x,y),{x' ,y')eQ. 

As for the bound on dniptJ '^^ have 9„ = cos 09a; + sin6'9j,, so that 

(3.8) dniptiani^) = -sin6'(s)a;'£;(y(s))(TXj.)(xB(y(s)),y(s)) 

rXE{y{s)) 

(3.9) -sm9{s) dy{TXu){x,y{s))dx 

Jx(s) 

(3.10) +cos0(s)(Tx.)(x(s),y(s)). 

Using the definition of x^, notice that if (xEiy), y) G Suppxi^, then 

either mi — > ;r or mt sup r^ > 1. 

ie/+ dy 2 iei+ ± dl,^±. 
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> 1. 



According to (3.6)-(3.7), the above condition becomes 

.,, . r |y-yik 1 ■ r {xE{y)-x^)± 

either mt — > - or mt sup ■; -. -. — —, — — r^ 

iei+ dy 2 iei+ ± \xE[yiTsm0{yi)dy) - Xi 

Since xe is a monotonous function near y — jji, this amounts to 

inf \y~yi\ > 6y/2. 

We infer that 

(3.11) 

The formulas in Appendix A imply that for all i G /_|_, in the vicinity of yi 

C 



{xE{y),y) eSuppxi' => V?: e /+, |y-yi|>y. 



(3.12) 
so that 



WEiy)\ < 



\y - y^\\^n\y - y,\\ 

{xE{y),y) e Suppx^ => Wsiy)] < CM{y) < 



c 



Sy\\nSy\ 



Recalling eventually that for all c > 1, for v sufficiently small, 

SnppdyXi. C IJ [xi- c5l._,x.i + cJ^ +] x [yi ~ 6y,yi + 6y] , 



ieu 



with Sl^j. = 0{\\n6y\ ^), this leads to the estimate on 9„Vt'|ao- We have indeed 

(3.13) ' 

XE{y{s)) r^Eiyis)) 



nXE(y[.S)) 

/ dyiTX^)ix,y{s)) dx\ < ^J \\dyX 



i^l|oo-'-|j^-j^i|<(5„ 



is) 



^\x — Xi\<c5^^X 



r^E{y{s)) 
+ zJ / X^\dyT\ix,y{s))dx 

( 4 

The estimate on 9s V'? is obtained by differentiating the identity 

rxE(y(s)) 

Jx(s) 



{TXi^){x,y{s)) dx 



with respect to s. Using identity (A.l), we infer that 

dsi}°\dn{s) = cos e{s)x'E{y{s)){Tx^){xE{y{s)),y{s)) 

+ s\ne{s){TXv){x{s),y{s)) 

t-XE(y(s)) 
- cos 6l(s) / dy{TXu){x,y{s)) dx. 

Jx(s) 



The estimates (3.12) and (3.131 above yield the desired inequality. 

The other trace estimates are derived in a similar fashion. Notice that is is 
hard to derive a sharp global estimate, since the angle 9{s) for s G Tw is in general 
different from the angle 9{s'), where s' G Te is such that y{s) = y{s'). Therefore, 
there is no simplification for terms of the type cos^(0(s))x'^(y(s)) when s G Tw 
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On the East coast, however, we can use the formulas in Appendix A, from which 
we deduce that 



cos'i9is))x'i,{y{s)) = 



cos 9{s) 



Vs G Te- 



Notice also that i/j^ vanishes on F^ by definition, that dnXp'^^j^^ is supported in 
{|i/(s) — yi\ > Sy/2}, and that all terms of the type 



t\rE 



XE(yis)) 

x{s) 



dyiTXu){x,y{s))dx 



are zero for s e Te- The upper-bounds for 9^9„f/;°|p , 9^9„V°|r follow. 

As for the jump estimates, we recall that with the notations of Chapter [2] 
Section |Z4l 



m 



2^B(yi) 



{TX>^){x,yi)dx, 



and theref ore t he jump of ■0? is of order one. Notice that the assumptions of Chapter 
[2] Sectic 



2.4 



imply that xi — x^{yi). The jump of the y derivative is given by 

[dy^th = -XE'iyi){-rxu)ix^iyi),yi), 
+xe {yi)iTx^){xE{yi),yi) 

dy{TXi,)ix,yi)dx. 



The points {x^{yi),yi) do not depend on ly. Therefore, either (a;^(yi),yi) is an 
East corner (1 e /+), and then {TX}y){x%{yi) , yi) = for all z^ > 0, or {x%{yi),yi) is 
not a corner, and then Xe (yi) is bounded by a constant independent of v. Hence 
the first two terms of the right-hand side above are bounded as z^ — ?> 0. Using 



inequality (3.13), we infer that 



dy{TXu){x,yi)dx 



< 



C 



5y\hi5y\ 



D 



At this stage, we have built on each subdomain Vt^ an interior term -0^ 

• which belongs to H^ (typically with derivatives of order 0{5^'^) in the 
vicinity of East corners), 

• and which approximately satisfies the transport equation dxijj'^ = t. 

Nevertheless "0? does not fit the boundary conditions, and therefore boundary layer 
correctors must be defined. Following the direction of propagation of the equation 
at main order, we start with the East boundary layers. 
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3.2. Lifting the East boundary conditions 

In this section, wc focus on an interval (si,Si+i) C T^- We recall that the 
domain [s]' , s^j^] C [si, Si+i] of validity of the East boundary layer is given by 

st := sup |. G (^.., '!l±^^ , ^i/3|0'| I cos^r^/^ > l| , 

^-+1 = inf |.s e (^fl^^, .,+i) , .i/3|0.| I cos^l-Va > i| . 

Easy computations based on the explicit rate of cancellation of cos 9 near s, provide 
then the following Lemma, which shows in particular that, because of the trunca- 
tion, the trace of dnipt is zero outside the domain of validity of the East boundary 
layer. 

Lemma 3.2.1. Under assumption (H2i), 

s+ - Si - C^v^^ . 

Under assumption (H2ii), 

, 4a / Inllni^l ^ / 1 

In i^ \ mv \mv 

Proof. Lemma [3.2. 1| is obtained by straightforward computations : 

• If (II2i) is satisfied, then for s > Si 

i^i/3|0'(s)||cos0(s)r^/3^CVi/3|s_s^|-^¥^ 
so that there exists a positive constant Ci such that 
(3.14) s+-s,^C,iy^K 

• If (H2ii) is satisfied, then for s > Si 

z.i/3|0'(s)i |cos0(s)r^/3 ~ cv/^!^^^!^^^ 

(s~ Siy 

so that 

-21n(.+ -.0+ ,, +" =-lnC'-^+o(l). 

-JWi - Si) o 

Computing an asymptotic development of sj' — si leads to 

, ia f In I In 1/1 ^ / 1 



D 



Using the formulas in Appendix A, wc infer in particular that in case (H2i) 

y{si) - y^ -^ Ciy^^ , 
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and {n + l)/(4n + 3) > 1/4, so that \y{s^) — 2/i| <C Sy. In a similar way, if (H2ii) is 
satisfied, 

\y{st)^y,\ ^ qinH-^exp(^a'^(^l + 6^^^+0(|lnH-^) 

~ q In z^r 1/2^,1/4 <<j^_ 

Therefore, in all cases, we have \y{sf) — j/i| ^ Sy, so that the extremities of the 
domain of validity of the East boundary layer are in the truncation zone. 

3.2.1. Traces of the East boundary layers. 

By definition of ^l'^, we have 

(^?)|r. = 0. 
We thus define the East boundary layer to lift the trace of 9„V'°, i.e. 

where 

\e{s) 



a\ 1/3 
COS B ^ 



V 



But then V'? + i^s'^i^T ^ez) does not satisfy anymore the zero trace condition 
on r^;. More precisely, since dnip^^Y vanishes for \y{s) — yi\ < Sy/2, we have the 
following trace estimate: 

Lemma 3.2.2. The trace ofipE satisfies the following bound 
(3.16) \Mz=o{s)\ < C- '^-——M{y{s)). 

\COSU(S)\'-''^ 

Moreover, its derivatives with respect to s satisfy 

\dsiJE\zMs)\ < C,y'/'S;'\cose\'/'M{y{s)), 

\d',^Eiz=,{-s)\ < C.^' W ^^'lfL,, M{y{s)). 

" |cosw(s)p/3 



These estimates are a straightforward consequence of Lemma 3.1.2 and of the 
formula defining -0b. 

Hence the remaining trace on F^ is non-zero, and must be corrected. Note in 
particular that the bounds above are too singular in the vicinity of Sj , i G /+ in 
order that we can lift the remaining trace by a simple macroscopic corrector. We 
lift the trace in two different ways depending on the value of s: 

• If s is far from any point Si,i £ /+, we lift the trace thanks to a macro- 
scopic corrector V's"^'^? which we construct in the next paragraph; 

• If s is in a neighbourhood of size one of any Si, we lift the trace of TpE 
thanks to the North and South boundary layer terms; we will explain the 
latter construction in the next section. 
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3.2.2. Definition of the East corrector. 

The idea is therefore to split each East boundary component [s,;, s^+i] in three 
subdomains (independent of v) 

• one which is far from the singularities 



'i+i- 



C\si,Si+i[ 



so that we have uniformly small bounds on the trace of ^p^^ on [a^ , cr~_-^]; 
• two which are (macroscopic) neighbourhoods of the singularities and such 
that 

(3.17) Vse [s„a+]U[a-+i,s,+i], |sin0(s)|>-, 

so that we can extend the North or South boundary layers on these parts 
of the boundary. 

We therefore define suitable truncation functions ip^ , (p^+i G C°°(R, [0, 1]) 

^tis) = 1 if s > a+, ^+(s) = if s < ;^(,s, + a+), 

(3.18) _ _ _ ^ ^ _ 

V>^+l{s) = 1 if S < cr,,+i, <P.,+i{s) = if S > l^iSt+1 + cr,+i). 



The East corrector is then expected to lift V'i ¥'i+iV'-E|z=Oi which is the part of 
the trace which is known to be small. More precisely, we set, for s G [sj_i, s^], 

(3.19) V^r''(s, z) = -V^s|z=o(5)^+(s)^r_^(s)xo(z). 
Using Lemma [3. 2. 2[ we infer that 

(3.20) ||V;--||^™(o) = 0(^i/3) Vm<3. 

Moreover, '/'^'i^qx p — 0, and dnip'^fg^ = by definition of xo- Notice also that 

(^0 + ^e{s, Xez) + ^r'")|[.^.r_ J = 0, 
9„(^?+^i=;(s,A£z)+^r'')|r. =0. 

Therefore, at this point, we have restored both boundary conditions on the 
subdomain [a^,a^_^] of [si,Si-|_i], but not the condition on the trace in the neigh- 
borhoods of the singularities. This is handled by the North and South boundary 
layers, which we now address. 

In the next section, we set 

(3.21) i, := ij^ + T/^l"-'- + ^e{s, Xe(s)z), 
so that 

dni'idnXFE = dnipt\dn\rE^ 
^|(s.,s.+i) = (1 - ftv7+i)i^E\z=ois), 

where (si,Si+i) is a connected component of F^;. 
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3.3. North and South boundary layers 

In order to lift the boundary conditions both on the horizontal parts and on 
the East/West boundaries close to the points {xi, yi) for 1 < i < fc, we then define 
North and South boundary layers. Without loss of generality, we focus on the case 
of South boundaries (the case of North boundaries can be deduced by a simple 
symmetry) . 

Denote by Si and Si_i the curvilinear abscissa of the endpoints of the horizontal 
part to be considered, with Si = Si-i if it is just an isolated point of the boundary 
with horizontal tangent. 

There are several point to be discussed in the present section: 

• In the first paragraph, we give a precise definition of the interval on which 
energy is injected in the South boundary layer (i.e. we define the support 



of the functions ^o and *i appearing in (2.11)); 

• The second paragraph is devoted to the derivation of regularity and mo- 
ment estimates on ips] 

• Eventually, we explain how we truncate "05 beyond the support of 'So f^nd 

3.3.1. Definition of the initial boundary value problem. 

We first define the extremal points of the interval on which some energy 
is injected in the boundary layer : 

• If cos 6'(s) > for s > Si (East boundary), we need to lift the East 
boundary condition on a macroscopic neighborhood of s^, namely [si, a^] 



where tr^ is defined by (3.17) 



Similarly, if cos0(s) > for s < Si_i, we define crj_-^ by (3.17). Truncation 
functions (pf,(p^_-^ are then defined by (|3.18 1. 



If cos 6'(s) < for s > Si (West boundary), we only need to lift the West 
boundary condition on a small neighbourhood of Si, precisely when the 
West boundary layer become too singular (Recall that the trace of ip^ is 
not zero in general in the vicinity of s^ on West boundaries). We therefore 
denote by a^ the arc-length beyond which no energy is injected in the 
boundary layer, and we expect that a^ > sf > Si, and \af' — Si\ <C 1. 

Note that, in order that the transition between the two types of 
boundary layers is not too singular, we have to choose \(jJ' — sJ\ as large 



as possible. But in order that the transport coefficient fe in (2.11 ) remains 
bounded in L^, we need that cos6'(s) is small for s < cr^ . We thus choose 
a^ satisfying the condition 

(3.22) / I cos 6'(s) I ds==iyi/^ 

J Si 

We further define (p+ e C°°(R, [0, 1]) such that 

(3.23) ip+{s) = if s < s+ and Supp(l - ip+) C (-oo, a+[. 
Similarly, if cos6{s) < for s < Si_i, we define a'~_^ by 



f ' '|cos6l(s)|ds = ^yl/^ 
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and (y37_i e C°°(R, [0, 1]) such that 

Supp(l - ^l_^) C]cr,ili,oo) and (/?ili(s) == if s > s^^^. 

For technical reasons, we also require that if (si,Si+i) C Tw and if 
cos 9 vanishes algebraically in a neighbourhood on the right of Si , 

(3.24) Suppa,^+ C (s+,s+ + h_^ 



(3.25) 



of course a similar assumption holds for f^_i. 
We further impose that for all fc G N, 



|9,V?IU - 0{\sf - afl"^) if (s„5, ± 1) c r 



E- 



Remark 3.3.1. There are two cases of South boundary layer terms which will 
not be considered in this section, since they both give rise to interface boundary 
layer terms: 

• when cos 9{s) < for s < Si_i (West boundary on the right) and cos 6{s) > 
for s > Si (East boundary on the left); 

• when cos 0{s) > for s < Si_i and for s > Si (East boundary on both 
sides), with Si_i < s^; 

We will indeed consider the corresponding South boundary together with the 
singular interface S in the next section. 



As in Lemma 



3.2.1 



we can compute asymptotic developments for aj' and (Ti_i. 

Lemma 3.3.2. The transition functions ip^ have controlled variations : 

• If cos > for s G {si,Si^i) (East coast), then there exists a constant 
c > such that 

and therefore 

ll9,V+lloo-o(i) yk. 

• IfcosO < for s e (si,Si+i) (West coast), then 

CT+ - s, - C^i/i^ 

if (H2i) is satisfied in a neighbourhood on the right of Si, and 

, 4a / 81n|lnj/| ^ .,, , ■ 
a+ - s, ^ 1 -^ L +0{\\ny\ 

if (H2ii) is satisfied in a neighbourhood on the right of Si. 
In this case, 

||9.V+lloo < Cu~^". 

Proof. The proof follows from the definition of af and ip'l on East and West 
coasts. We check in particular (see Lemma 3.2.1 1 that we always have a'l > s^ . In 
the case of a West coast, using assumptions (3.251, we have 

lia^V+lloo = o(.-^-/(^"+-^)) (resp. iia^V+lloo - o (il^)) 
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if assumption (H2i) (resp. (H2ii)) is satisfied in a neiglibourliood on the riglit of s^. 
Notice that the most singular case corresponds to (H2i) with n = f , which leads to 
the estimate of the Lemma. D 



s> a*: i|j^ compensates 
the traces of LJJ, 5 ijj 




"Mixed" zone 



South coast: 
s< s*: ip compensates ''O^ ^"0, sin 9=-1 

the traces of Tjj, 5 ijj 



Figure 3.2. The interplay between the North/South and 
East/West boundary layers 



Define 

(3.26) 7.-i,.-(l-^+)(l-¥'r-i)- 

The South boundary layer is therefore described by 

ds^s + h(s)dz^s + /i(s)9|V's = 0, s> a- 1, Z > 0, 

V'Z|X=0 = *0, dz^-'s\Z=Q = ^1 

where b{s) ~ — i^^^/^cos0(s)/sin0(s), ^{s) = — l/sin6'(s), with boundary condi- 
tion 

'^a{s) ■■^-li-i,S\dn{s), 
(3.28) , 



and zero initial data prescribed at s = (7^_i- We recall that ?/) is defined by (3.21 1. 
Note that, when the two types of boundary layers meet (that is on the support 
of (7i-i,i)'): the width of the North/South boundary layer is much larger than 
the one of the East/West boundary layer. Indeed, the width of the North/South 
boundary layer is always v^'^ , while that of the East/West boundary layer is, using 
Lemma |3.2.1[ 

/ \ 1/3 ( 1+1 

-1 (,+\-^^J^_\ <r} ^'"+' incase(H2i), 



^E,w^^i) - \ ^^^,air.+ \] -^1 ,,i/4n„ ,,1-3/2 



z^i/4|inzy|-3/2 in case (H2ii). 
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Therefore, it seems more accurate to talk about a superposition of the boundary 
layers, rather than a connection. 

3.3.2. Estimates for 'ipN,s- 

Lemma 3.3.3 (Trace estimates). The functions defined by 

^o(s) := -7j-i,i'(A|an(s), 

*l(s) := ~l^^^*J^-l,^^n'^lJl9n{s), 

for s € (cr,^]^, Si+i) are such that 

||a,"*o||L- + ||a,^^*i|lL~ < Ciy-""/' forke{0,--- , 3}, 

\\ds-9o\\L^ + \\ds-9i\\L^<C, 



\d',^„hr + m^i\\L^<cd-\ 



la^vi/ollii + 1193*111^1 <C6;',^-'/^. 



Moreover, 



||6*oIUi + II^'*i||li=0(1), 

Proof. We recah that ij; = Tp° + ■0^°'"'' + V'_e(s, Ab(s)z). Notice that thanks to 
assumption (H3), tpt i^ continuous (and even C^) on the interval {iJ~_-f^,a^). The 
estimates are slightly different depending on whether (si_2,Si_i) and (si,Si+i) are 
portions of Te or Tyy- We focus for instance on the portion {a^_-^^, Si^i), keeping 
in mind that the portion (sj,cr^) is analogous. 
Connection with West boundaries. 
If {(J^_i, Si-i) C T\y, then by definition of (t^_i, 

(3.29) Vs e ((Tr„i, s,_i), \yis) - y,\ < v^'^ « by. 

Moreover, ?/! = ^/i^ on {oJ_^^ s^) in this case. Therefore Lemma ', 
of the trace estimates on ip^ stated in Lemma 



3.1.2 



3.3.3 



IS a consequence 



First, the L°° bounds together with the definition of ip^_^ huply that ||*oI|l=° = 
0(1). Moreover, for s e {(j^_j^,Si), using Lemma 3.1.2 and inequality (3.291, we 

have 

|cos6i(s)r 

Sy\\nSy\ 



\dsrt\ani^)\ <C(l 



so that, using the definition (3.17) of af together with the definition (3.51 of Sy, 

rt 

II 5«V'?|ao II Li(.r,, so < C(l + z.-i/4|lnH-^/^(ln|lnH)M \cos9{s)\ds 

< C(l + |lni/|-^/5(ln|lni^|)'3). 
We infer that 

\\%\\L^i.-_^.,s.) < Ut\an\\oo\\ds^:h^ia^_^,sO + \\ds^t\on\\LH.r_^.s,) < C 
In a similar fashion. Lemma 3.3.2 and assumption (H2) yield 

< ^^/^||a„^?|,ollL~(.,:^^,..) 

1 



l^ll 



L°'(crr 



< Cv''H\ + 



5y\\n5y\ 



<c, 
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and eventually 

1 I cos 6*1 



/ j,l/4 j,l/2 X 

The higher order estimates are obtained in a similar way. We use the estimates 
of Lemma |3.3.2[ namely 

In a similar fashion, 

implies that for s G (Cj^x' ■*«) 

^.2 , 



Lemma 



3.1.2 



Wiis)\ + Wiis)\ < ^(^;^(l + ^)+--^/^(l + ^^)) 



|5,^vi>o(,)| + |a,^vi/i(s)| < C5-M1 



y / \ y 

+c(|5,vi + l9.Vr-il), 

cos 01 I cos 01' 






+C.--" (l + ^^^) + C {\dU\ + |5.V.^rl) 



from which we easily infer the estimates of the Lemma. 

The estimates on 5^o and 6^1 are a consequence of the following estimates on 

b: 

6 = on (sj_i,Si), 

by definition of CTj^j^ on West coasts, and 

Indeed 

Using assumption (H2) together with the definition of cr^_i, we infer the desired 
result; notice that the most singular case corresponds to rt = 1 in (II2i) for the 
estimates on b and dgb, and to n = 2 for the estimate on d^b. 

Connection with East boundaries 

If (cr,^-^, Si_i) C Te, the estimates are different in several regards: 
• The function 7i_i,i has bounded derivatives; 
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• Because of the truncation Xv^ the traces of tp^ and dni^^ are identically 
zero on the vicinity of Si_i; 

• The trace of V'? is zero on [<7~_-^,Si-i). 

• The normal derivative of "0 = ■0^ + tpE + 4''^^^ is identically zero on 
(cTj^^, Si_i) by definition of ^e, so that ^i = on [cr^^n Si-i)] 

Therefore it suffices to prove that ^i-i,i4'E satisfies the desired estimates on the 
interval {iy~_ij Si-i). Notice that the estimates on (sj_i, s^) can be treated with the 
same arguments as in the first case. 

We use the estimates of Lemma 3.2.2 and assumption (H2), together with the 
formulas of Appendix A, in order to compute y{s), cos0 and 6' in terms of s. We 
infer that on (a~_^, Si-i), ^o satisfies the following bounds: 

• If (H2i) is satisfied in a neighbourhood on the left of Si_i, then for k — 

1,2,3, S £ (cr^^lpSi-i), 

l*o| < C^'/' {\s - s.-il-^l ln(s,_i - s)r'lc5y(-.+i)<|,_,^_^|<i/2 + l) , 
We infer in particular that 

Since „1"+'?, < 4/3 for all n > 1, we infer that 

3(n+l) ' — ' 

As for the other estimates, we have, for k — 1,2,3, 

It is easily checked that the most singular estimates on the L°° norm, and 
on the L^ norm as soon as fc > 2, correspond to n = 1. We obtain 

\\ds^"hH.r_^,s.^.)<C\\n,.\-\ 

Notice that these estimates are always of lower order than those derived 
in the first case. 

• If (H2ii) is satisfied in a neighbourhood of the left of Si-i, then for s G 
((Tj^i,Si_i) and k = 1,2,3, 

|*o(s)| < Cj^^^^\s - s.^il'^exp ( _"^_ I j hy(^)-y,^,\>Sy/2, 
\d^Ms)\ < C.'/%'\s - .._i|-^-^'= exp ( 3 "^_ ) l|,(.)-,._,|>.„/2. 
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Moreover, in case (H2ii), we recall that 

y{s) - Hi^i -- C{s- Si„i)^exp ( -■ 
so that for fc = 0, 1, 2 



s- Sj-l| 



and similarly 

R'*°l!L~K:i„.._O^C'|lnH^^/^ 

Therefore we retrieve the desired on the norm of $o- Concerning the 
estimates on &^0i we use the following inequalities on (Tj^x' ■^j-i) 

15.(6*0(5))! < C,^'^'^(^5y'\cos9\'/' + \e'{s)\\cose\-'/^)Miyis)) 
mbMm < C.Vi2(5-i|cos0p/3|0'(.)| + ^rai_)>i(y(,)). 

Once again, we distinguish between (H2i) and (H2ii), and we obtain 

||&*o||li <C\ln,^\-\ 
||3.(&*o)||li < Ciy^^^^d;^ < Cv-^/'^\lnu\-'/'^\ln\lnv\f, 

D 



By Lemma 2.3.1 we then have the well-posedness of equation (3.271, as well 



as suitable a priori estimates for the solution i/is 

IIV'5(s)|U^(R^)<C(|vI.o(s)| + 1*1(5)1) 

+c r (Ms')+i6(5')i)(i*o(5')i+i^i(5')i)rfs'- 



In order to show that the South boundary layer is stable, we will need additional 
estimates on the boundary layer terms, namely regularity and moment estimates. 

Lemma 3.3.4 (Regularity and moment estimates). Let ^s ^e ^he solution to 



(3.21) with boundary conditions (3.28), defined on 

{s>a-i/V.s'e[ar_^,s],sin0(s)^O}. 
Then 

11(1 + Z'=)V's(s,Z)Lo.(i|) + 11(1 + Z'^)dl^s\\Li^ = 0(1), 

(3.31) 11(1 + Z'^)ds^s\\Lr(Ll) + 11(1 + Z^)dlds^Ps\\Ll, = 0{5-^), 

11(1 + Z^)dl^Js\\Lr(Ll) + 11(1 + Z^)dldl^s\\Ll , - 0{u-^'^\ InH-^). 
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Proof. Such a priori estimates are obtained by combining precised energy 
estimates for the parabohc equation (3.38), together with high order estimates on 
the traces. 

• Define as in the previous chapter 

S{s, Z) := -*'o(s)(Z + 1) exp(-Z) - *'i(s)Zexp(-Z) 

+ 6(s)*o(s)^cxp(-Z) - 6(s)*i(s)(l - Z) cxp(-Z) 
- ^(s)*o(s)(^ - 3) exp(-Z) - ^i{s)■^^(s){Z - 4) cxp(-Z). 
From the trace estimates in Lemma [3. 3. 3[ we deduce that, for all fc > 0, 

II^'^IIli<c, 
Wz'^OsSWl^ < C5y\ Wz'^dlsWLi < c5-\-^i\ 

Notice also that by definition (3.18), ( 3.23| ) of kp^ , we have Supp(l — '-Pl_^ C 
\(7~_^-, oo[, SO that ^0 and ^I^i are identically zero on a neighourhood of (y^_i- There- 
fore 5,Vs|.=^r_^ = 0. 

• In order to derive estimates on 9^'0s for fc > 1, we lift the boundary conditions 
to get 

dsQ + Ks)dzg + li{s)d%g = 5(s, Z), s e (0, T), Z > 0, 

(3.32) 5|,=o = 0, 

g\z=Q = 0, dzg\z=o = 0, 

then differentiate the equation and proceed by induction on fc. We start with 

dMa) + d^{h{s)dzg) + ^';{^i{s)^%g) = d^S . 

Integrating against 9^g, we get 






3=0 ""-i 

and therefore, since /i is bounded from below, 

+cY^\\d':-^b\\,.,.^^-^^^^\\dig\\'^l^^^^^^^^^^ 

3=0 

fe-1 



3.3. NORTH AND SOUTH BOUNDARY LAYERS 57 

Eventually, using the Cauchy-Schwarz inequality, we obtain for k < 2 

5fc„l|2 I II afcQ2 „||2 



(3.33) m9\\U^^^-_^,^^,.^, + \\d^dh\\Ui^-_^,hLl) 

|2 



(3.34)< C||a,^-5||2 , 



fc-i 

/ ^ I « ^iiL°°((T._j.s) II * ^^iiL^(((T._j.s)xR+) 

Next recall that ^ is bounded, and that all its derivatives are of order 0(1) 
/^ > Mo > 0, \\difi\U < C. 

In order to bound the term involving b, we have to distinguish between the cases 
when (si_2,Si_i) C Te or Tw'- 

• If (s.i_2,Si-i) C Tw, then fe = v^^/^cosB^ji satisfies 

Notice that in this case, we also have (si,Si+i) C Tw since we have 
excluded for the moment the case leading to a discontinuity line S. 



Inequality (3.331 therefore yields, for fc = 0, 

\\9\\lt(lI) + II<9|5I1l1, < CII^IIl^lI) < c, 
then, for fc = 1, 

\\ds9\\Lr(Ll) + \\dldsg\\Li^<C5y\ 
and, for k ~ 2, 

If (si-2,Si-i) C r^, then for s G (crj^j,Si_i) 

for all j and for all s > Cj^]^. In this case, we refer to the proof of Lemma 
13.3.31 and we state that 

\\Shm.-_,,s.^,),Ll)=0{\ln,.\-'), 
Inequality ( 3.33[ ) therefore yields, for fc = 0, 



then, for fc* = 1, 
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and, for k = 2, 

There remains to derive bounds for s > Si_i. There are essentiaUy two 
cases: 

— either (si,Si+i) C T\y, and in this case the estimates on (si_i,cr^) 
are similar to the ones of the case (si-2, Si-i) C Tw- The only change 
lies in the initial data: we merely replace 

by 

P'5|s=s,_J1l| =0(z/-'=/4|lnz.|-i) forfc = 0,l,2. 
We infer that 

WahriLi) + II.9IIl2_^ =0(1), 

\\dsg\\LriLl) + \\ds9\\Ll,^0{S;'), 

Wd^ghriLi) + llCgLj,, = 0K^/^| in^-'). 

— or {si, Si+i) C r^;. Notice that in this case, Si = Si-i (i.e. there is no 
flat horizontal boundary), otherwise the corresponding piece of flat 
boundary would be the beginning of a E layer, which will be treated 
in the next section. Whence we can treat the interval (s^, a^) exactly 
as {a~_i,Si-i), and we obtain 

II.9|Il»(l|) + II.9|IlJ,, =0(|lnt.ri), 
\\d'.9\\Lr(L%) + \\d^9\\Ll, - 0{i.-'^^\ lnH-^)for fc = 1,2. 
• We get moment estimates in a similar fashion 

dsiZ'^g) + Z\h{s)dz9) + Z\ii{s)d%g) = Z'^S. 
Integrating against Z^g^ we get 
1, 



--\\z'g{s)\\i.^+ ^.{s')\\z'^y{s')\\i.^ds' 

< r \\Z^S{s')\\^.^\\\\Z^g{s')\\^.^ds' 

+ 2k f \h{s')\\\Z^g{s')\\L.^\\Z^-'g{s')\\^.Js' 

+ 2fc r |A.(5')|||Z'^a|g(s')L|ll^'-'azg(s')||L|rfs' 
+ k{k-l) r |Ms')lll^'a|g(s')Ld|Z'=-2g(s')|L2ds', 



A simple recursion then leads to the expected uniform bound for the moments of 
ips and dlips- 

• Combining both arguments we obtain the moment estimates for higher order 
derivatives. 

D 
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3.3.3. Extinction and truncation. 

We emphasize that ''Ps\s=a+ ^^ ^"-"^ smaU in general: there is still some energy 
in the South boundary layer, even though the boundary conditions are absorbed by 
the West boundary layer only for s > aj . We deal with this somewhat unexpected 
phenomenon by propagating the boundary layer term ips beyond af . We prove 
that when cos 9 becomes bounded away from zero, ijis vanishes. As before, there 
are two cases, depending on whether (s^, s^+i) is a part of Tw or Te'- 

• First case: (si, Si+i) C r^: looking closely at the proof of Lemma 3.3.4 we 
see that in this case, 

(3.35) Us\\L^ii,sa),Ll)=0{\\nv\'^) 

for all fixed s > af (independent of v) such that | su\9\ > 1/2 on [si, s). Therefore 
ips is small on the West boundary layer, on intervals of size one. 

• Second case: [si, Si+i) C T\y'- we use the following result: 

Lem ma 3.3.5 (Extinction of th e Sou th boundary layer). Let ■05 he the solution 
to {3.21) with boundary conditions (3.28), defined for s > (t^_i as long assn\6{s) < 
0. Then 



.-1/4| 



',e{s')\Us{s')\\l.ds'<c. 



(3.36) 



Before proving Lemnia [3.3.5[ we now explain how we truncate the South bound 
ary layer on (s^, Si+i). We introduce a function 7i-i.i G C°°{dVt) by 

Supp7i_i,i c]si_2,Si+i[n{sin6' < 0}, 

7i_i i = 1 on a neighbourhood of Supp(l — i^f) n Supp(l — '■P^-i)- 

Notice than we can always assume that 7i_i,i has bounded derivatives. 
The South boundary layer on (si-2, s^+i) is then defined by 

li~i.i{s)^s{s,zv^^/^). 

Using either Lemma 



3.3.5 



or estimation (3.35), we infer that II V's II l2 ^Odlni^l '^) 
on the support of T^TTT^ We will prove in the next chapter that this is enough 
to ensure that the error terms generated by 7i-i,i are admissible in the sense of 
Definition [T321 



Proof of Lemma [3.3.51 By Lemma 3.3.4 for every s G (si,Si+i) such that 
sin 6'(s) < 0, 

sup \\{1 + Z)^s\\l-=0{1), 



sGFi-i'^l 



and 



(l*o(s)l + WM\ + (a^(«) + IH*)l)(l*o(s)| + 1*1(^)1)) ds = 0(1). 



Moreover, for s > af , ^o(s) = ^i(s) — 0. Litegrating (3.27) against ZipSi wc infer 
that there exists a constant C > 0, which does not depend on u, such that 

Zi}l{s)dZ+ [ iy-^/^\ta.n9{s')\-^\\^Psis')\\l2ds' <C. 



Since | tan 6*1 ^ > | cos 6*1, we obtain the desired estimate. 



D 
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3.4. The interface layer 

Let us now focus on the interface S = 90+0357^. We recall that, by construc- 
tion, V'? is discontinuous across S so that 

• we need to lift this discontinuity to get a smooth (i?^) approximate solu- 
tion, 

• we then have to define a corrector (in the form of a singular layer term) 
to remove the energy created by the lifting term. 

3.4.1. The hfting term ip^'^K 

We first construct the term f/)''^* which lifts the discontinuity across E. The 
construction has already been explained in Chapter vA Since a^ is now defined 
precisely by (3.22) for West boundaries and ( 3.17[ ) for East boundaries, we can 
define a^ in the following way 



(3.37) 



So 



if/i = {si}, 

if /i = [so,si] and (cr(7,so) C Tw, 

if ^1 = [so,Si] and {crQ,So) C Te- 



Notice that if /i = {si}, then we always have (crj^,si) C Tw (an isolated point 
of cancellation of cos 6* on the East boundary does not create a discontinuity line). 
The reason why we choose to take sq instead of cfq as an initial point in the last line 
is essentially technical: we have chosen to write the equation on the discontinuity 
type boundary layer in cartesian coordinates, rather than curvilinear ones (see 



equation (2.191). This is legitimate as long as cosO <^ 1, so that the normal vector 
coincides almost with —Cy, while the tangent vector coincides with —e^- However, 
by definition of Uj on East boundaries, cos6 takes "large" values on {a^jSo), 
and therefore the present construction cannot be used. We by-pass this problem by 
considering separately the East boundary for s < sq and the discontinuity boundary 



layer (which will be justified in Remark 3.4.4 1 
We thus define a truncation ip by 



'P 



if/i={si}, 

if /i = [so,si] and {(Tg.sa) C Tw, 

if /i = [so,si] and (ctq^sq) C Te- 

Notice that in the last case, '0? vanishes identically in a neighbourhood of (xq, yo), 
and thus there is no need for a truncation in s. 



We then take (see (2.16)) 



V'"''(^,2/) = X+X ( ^^ ) Hx) + b{x){y - V,)] 



where x is a truncation function as in (2.17) and a and h are defined by 



b{x) = -(l-(^(s))sin6l(s)a„V't%o(s) 

+ (1 - Lp{s)) co&e{s)dsi}^i\Qi^{s) - co&e{s)ip'{s)'ijj°^Q^{s), 

a{x) = -(l-(^(s))V'°|an(s)-^(2;(s))(y(s)-yi)' 
for Xi < X — x{s) < x{a^), and by 



b{x) = -[a,v?]s 



for X < x^. 
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For X > x{a^), we merely take a{x) == b{x) = 0. 

Eventually, in order to simplify the present analysis, we truncate the function 
T in the vicinity of inf /i and sup /i , even in the case when these points do not 
belong to /+. Of course, this does not change anything to the estimates of Lemma 

inn 

3.4.2. The interior singular layer ip'^. In the previous chapter, we have 
described how the singular layer is connected to some North/South boundary layer 
at the East end. Assuming that the junction point {xi,yi) is some South-East 
corner, meaning that [so,si] corre spon ds to some South boundary and [si,S2] to 
some East boundary (as in Figure 2.2), the singular term is given by il;^{x,Y) = 
(t)^{x,Y -Y^{x)), where 



4 J,S 



(3.38) 



= -,5T"f*(x,y_(a;) + Z), a; < a;(cr-), Z > 0, 



^\x=x{a-) ~ ^' 



■'\z=o 



= d. 



z<Pz=o 



= 0. 



where (5t''^' is defined in Lemma 2.4.3 and Y- is defined by (2.20). 



The energy estimates for North and South boundary layers given in Lemma 



|2.3.1| can be extended to account both for the geometry and for additional regularity, 
which leads to 



Lemma 3.4.1. Let 0^ be the solution to (3.3S), and ip^{x,Y) = 0^(a;, F 



y_(x)). 

Denote by {xY tVi) — (2;(s]^), y(sj^)) the projection of {xi,yi) on the West 
boundary. 

Then, for all x < x^ , 



1^^ 



L°=((s,£c((T(7)),L2(y_^oo)) 



\d.A'' 



lL''((s.x(<T-)),L2(y_,oo)) 

Furthermore, for all integers k, I with fc + I < 3, 



\dl^^^ 



0(1), 

0(.-Vl9)^ 
0(^.-2/19). 



, , ■r-i l + 2I + 8fc 



Remark 3.4.2. The function ip^ is defined beyond the West end of {{x,y) G 
il) / y ~ ?-}•' indeed, for some values of s in a neighbourhood of s^ , x{s) may 
be smaller than x^ (notice that this is the case as soon as sin^(s]^) ^ 0, i.e. as 
soon as the tangent to dfl at s = s^ is not vertical). Thus the trace ip^ involves 
values of tp^ at points x < x\ 



w 



Proof. Such a priori estimates are obtained by combining precised energy 



estimates for the parabolic equation (3.38), together with a refined description of 



the geometry IC, and high order estimates on the source terms. 
• Estimates on (b^ . 
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Denoting hy S — S{x,Z) the right-hand side of ( |3.38 ), we recaU that 
5(.T, Z) = (1 - l,<,,,z+y_(x)<o) X{Z + Y^ix)) [a'(x) + iy'/^b'{x){Z + Y^{x)) 
~ {^-i.<.uZ+YAx)<o)/^Hz + Y^{x)) [aix) + iy'/Hix)iZ + Y^{x)) 
- 4 (1 - l,<,„z+y_(,)<o) X^^'XZ + Y^{x)y/^b{x). 

Notice that the indicator function does not create any singularity on the domain 
of integration x € {x,x{a^)) x (0, oo): indeed, x ^^nd x are identically zero in 
a neighbourhood of zero, so that ly<ox'''^''(^) i^ i^i f^-^t a C°° function for k > 2. 
Moreover, for x < xi, a and b are constant, so that 

lx<.il(z+y-(.))<o X{Z + Y^{x)) (^a'ix) + j.i/V(x)(Z + r_(a;))) 

Eventually, for x — xi, F_ — 0, so that lz>olz+y_(xi)<o = 0. 
It is easily checked that 

\\Six, Olk^tR,) < Ci\a\ + \a'\ + z.V4|5| + ^1/4|;,.|)^ 
SO that, using Lemma |B.1| in Appendix, 

II'S'||li((2,x(<j-)),l2(r^)) < c. 

According to Lemma 2.3. 1[ we then have 

110 ||L~(2,a;(cr-)),L2(R_|_)) + \\dz<l> \\l^({xMc-))>^R+) = 0{l). 



0. 



In order to derive estimates on dl^ip^ for A; > 1, we differentiate equation (3.38) 
and we proceed by induction on k. Notice that because of the various truncations 
and of the choice of a~ , (f^ is identically zero in a neighbourhood of x = x(cr~), 
so that the initial data is always zero, together with the boundary conditions. 
Moreover, as explained before, the indicator function in the right-hand side does 
not raise any singularity. Using the same kind of estimates as in Lemma |3.3.4[ we 
infer that 



(3.39) 






fe-1 



+C^\\Y_ 



IIl-/3(2,.(.-)) (||ai</.^Ljo(^|) + Wdldlcf' 



• Estimates on d^S. 



As in the proof of Lemma 



3.3.2 



the most singular estimates for Y_ ' are ob- 
tained when cos 6 vanishes algebraically near s = si, with the lowest possible ex- 
ponent n. Moreover, by assumption (H4)), the exponent n above satisfies n > 4. 
Hence we have 



(3.40) 



rilco<C|lni.|i/5i.-^-/20^ 



\Y^-^\W. 



< C|lnj/|i/5j/(^ 
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We now use the estimates of Lemma B.l in Appendix B together with (3.40) 
We have 

\\d^S\WL- < C{{Ai+A2) + \\YL\\oo{Ao + Ai)), 

\\dlS\\LiL'^ < C{A2 + As) + \\YL\\o.{Ai + A2) 



-C{\\Y!! 



\YL 



J(Ao + Ai), 



l!5^5|Uii2 < C{A3 + A4) + \\YL\\oo{A2 + A3) 
+Ci\\Y::\\^ + \\YU\l)iA,+A2) 



so that 
(3.41) 

• w^'°° estimates on ^^ 



YH 



|a,^5||^.(^i) < c-'/''. 



YL 



KAo + A,), 



• Gathering (3.39), (3.40) and (3.41), we prove easily by induction that for 

fce{0,---3}, 



\<-'z"x'r \\L^{(x,x(a-),L^{Yl+)) 



The first estimates on -0^ foUow from the identities 

dlil^'^{x,Y) = dl4P{x,Y-Y_)-2YLd^dzcjP{x,Y-Y_) 

-Y!!dz4>^{x,Y- y_) - Y!_^dl(t)^{x, Y~Y^). 
Then, notice that for x < x{t'^), Y^ = Y^, so that 

Wr llH"'-.°=((2,3;(t+))x(?_,oo)) = W^ lwfc-=°((S,a;(t+))x(0,oo))- 

We then use Sobolev inequahties in order to bound the L°° norms. More precisely, 
we use repeatedly the following Lemma, whose proof is postponed to the end of the 
paragraph: 

Lemma 3.4.3. There exists a constant C > such that for any f G _ff'*(R_|_) 
satisfying f\z=o — dzf\z=o ~ 0, we have f G W^'°° (R-^) and for fc = 0, 1, 2, 3, 



Notice that for x < x{t^), (p^ satisfies 

d^<jp - d%<jP = {iz+Y-<o - i)x^^\z + Y-){a + v''%{z + y-)), 

so that for k > 2, 

dl<f = d%dl-^<f, x<x{t+),Z>Y-. 



Lemma 3.4.3 and equation (3.381 imply that for k> 1, 



\d: 



ife^Si 



L°°((2,x(t+))x(0,oo)) 



< 


\\Ox<P \\l^(lI)\\"z0x9 IIl~(L|) 


< 




< 


7k fc + 1 
(Jll 8x19 J/ 8x19 



< C'ly 19 8X19, 
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while 

||q; j,S|, < ^11 ,E||(7-2/)/8||^4 ,E||(l+20/8 

\\Oz<P llL~((£,2;(t+))x(0.oo)) ^ ^WP \\l^{lI)\\"z<P IIl~(L|) 

(l+2Z)/8 



< c(||a,<^^L^.(i|) + ||5L^.(i|)) 



_ 1 + 21 
< CV 8X19. 



Similar estimates can be derived for terms of the type d^d^'^cj)^, with k + l < 3. This 



concludes the proof of Lemma 3.4.1 



D 



Remark 3.4.4. • Because of the function y_ and the coefficients a and b, 
the profile ip^ depends on v in general, which is not entirely satisfactory for an 



approximate solution. However, the estimates of Lemma 3.4-1 show that ip^ is 
uniformly bounded in L°°{L'^) n L?{H'^). 

Let X G {x,Xi) arbitrary. Then x{tf) > x for v small enough, so that Y_{x) = 
Y- for X e {x,x). Extending i/j by zero below Y-, we infer that ip is uniformly 
bounded in L'^ {{x , x) , H'^ (R.)) . Up to the extraction of a sub-sequence, ip'^ converges 
towards a function i/i^ weakly in L"^ {{x , x) , H'^ (R)) . Passing to the limit in the 
equation satisfied by V'^j we obtain 

d^ij^ - d^i;^ = S{Y), in V'{(x,x) x R), 

where 

r^tivi) 
S{Y)^~dx'^^\Y),witha = - / T{x,yi)dx. 

■'^EiVl) 

• Notice that once V''''* and "tp^ are constructed, there remains a non zero trace 
on the East part of the boundary (si, S2), and also on (sfc_i, Sq) if {s^-i, Sq) C T e 
(notice that a~ = Sq in this case according to (3.371, and therefore the trace of ipE 
on (cTg'jSo) is not lifted by ip ). 

These traces correspond to 4'E\z=o{^ ~ fi'P2) '^'^'^ ^"£12=0(1 ~ Va'Pk-i)- They 
can be lifted thanks to South type boundary layer terms, exactly as explained in the 
previous section. 



We conclude this paragraph with the proof of the Sobolev estimates used in 
Lemma 13.4.11 

Proof of Lemma 13.4.31 The estimates on / and /' are derived in a classical 
fashion: for instance, write 



/•OO 

/(Z)2 = -2 / f{Z')dzf{Z') dZ', 
J z 



I z 
so that 

||/||L<2||/||2||9z/||2. 

Integrating by parts and using the fact that /|z=o = c^z/|z=o = 0, we infer that 

\\dzf\\l^-l fdlf 
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and similarly 

(3.42) WdlfWl = I fd%.f. 

Using several times the Cauchy-Schwarz inequality leads to 

ll/lloo < ^/2||/||^/«lla|/li^/^ 

The estimate on dzf goes along the same lines. The estimate on d%f is a little 
more tricky: we write 



\dlf\L < 2\\difh\\dy\\ 



2 



l/2||o3 ^11 ||C)4 i.||l/2 



(3.43) < 2||/||^'1|ai/||2||5|/|| 

On the other hand, an integration by parts yields 



2 



and 



14/112 = - / dlfdy - a|/|z=oa|/|z=o 

JR+ 

jRi 



Therefore, using (3.421 



PI/1I2 < c\\f\\'/'\\dy\\i/' + c\\dy\u\dy\\i^'\\dy\\i/'. 

We deduce that 

Wdifwi < c\\f\\l^'\\di,f\\i/' + c\\dy\\ti'\\dy\\i/\ 

Inserting the above inequality into (|3.43 1 and using Young's inequality leads even- 



tually to 

\\dy\L<c\\f\\i/'\\dy\\i/\ 

so that 



d'zfh < currmn 



2 



The last estimate on ||9|/||oo follows easily. 

D 

3.4.3. Connection with the West boundary. The question is then to de- 
scribe the West junction, that is the connection with the West boundary layer. The 

traces of the discontinuity terms -0''^' -I- V'^X I ^~s^^ ) on the West boundary create 
yet another singularity in the West boundary layer term: indeed, the coefficient 



A{s) in (2.7) now lifts the traces of V"™* + "0 , and thus changes abruptly near the 



int „/,0 I ,/,litt 



point sf e Tvi/. We recall that V'"* = V't + V" 

At this stage we encounter another difficulty: the function dsA can be written 
as 



1/ 



1/3 



dsA±is) = /3± ^ ' cos0(s)sin0(,s)92^™*(a;(s),y(s)) 
cost^ " 

-|-differentiable terms. 
The function dy{'tpf + ?/'''^*) is discontinuous across E: more precisely, 

and the jump of dyijjf across S is constant along S and of order 6~^. As a conse- 
quence, A' is not continuous, which is problematic, because the error estimates for 
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the West boundary layer involve 17 norms of A" , as we will see in the next chapter 



(see the proof of Lemma 4.2.4) 



In order to avoid the apparition of a Dirac mass in A" , we add yet another 
corrector term which has the following form 

(3.44) r- . Amhly>yAy-V.fx (^) X {^^^^) • 

It is obvious that il)^" £ H'^{fl). Moreover, ?/)g""'' satisfies the following estimates: 
Lemma 3.4.5. For all fc € N, we have 

Moreover, 

Proof. The estimates are fairly straightforward, noticing that for I € {0, 1, 2}, 

fceN, 



0[v''^). 



Therefore, recalling that 5y = i/^/^| lni^|^/^| In | hiv\ \^^ for some /3 > 0, 
RVriloo=0(<5;V/2|lnH'=/') = 0(|lnH'^|ln|lnHI'^), 

The W^'°° estimates follow easily. Moreover, since the support of ij}'^'^^ is included 
in the rectangle 

we infer that 

ll9.^rnU^(o) -0(|lni.ri/5|ln|lnz.| p/5|lnz.ri/iOj.i/8) ^o(^i/8)^ 

D 
In the rest of the paper, we set 

By definition, '0s is localized in the vicinity of S and such that 

• Vt° + -02 belongs to H'^ift), 

• the energy created by V's on Vt^ is negligible. 

Note further that, by assumption on S, the trace of -02 + V'? + i'N,s on the 
boundary Tw vanishes in the vicinity the points Si € Tyy- 

Remark 3.4.6. If one looks carefully at the construction above, there are several 
reasons why the assumption (H4) is necessary: some of them are technical, and 
others are probably more fundamental. We first notice that if Ii = [so,Si] with 
Sq < Si, assumption (H^) is always satisfied, as dVt is C* near sq- Therefore we 
focus on the case Ii — {si}. 
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• First, assumption (H4) is required so that the coefficients a and b have 
derivatives of order four; if cos 6 vanishes only at order n nearsi, withn < 
3 in (H2i), then such a regularity is not achieved in general. Therefore, 
if (H4) is not satisfied, an additional regularization of a and b must be 
performed. 

• Moreover, the estimates on tp'^ are more singular if the cancellation of 
cos 9 near s — Si has a lower order. This could have two different conse- 
quences: on the one hand, it could be possible that the estimates on tp^ and 
its derivatives in the interior offl are not sufficient to prove that the error 
terms are admissible. This, however, is unlikely: indeed, the estimates are 
exactly the same as the ones for North and South boundary layer terms, 
for which we are able to prove that the error terms are admissible even 
when cos 9 vanishes at the lowest possible order. 

On the other hand, if cos 9 vanishes at a low order near s = si, the 
estimations on ||V'^||h"='°° become much more singular, and therefore the 
trace of ip^ and its derivatives on the West boundary are much larger. 
In this regard, assumption (H4) does not seem to be merely technical: it 
is possible that the discontinuity boundary layer may destabilize the West 
boundary layer if ip^ becomes very large. This could in fact be achieved 
with a low order of cancellation of cos 9 near si. However, such consider- 
ations go beyond the scope of this paper. 

3.5. Lifting the West boundary conditions 

We recall that the boundary layer term on West boundaries takes the form 
■0w(s, Z) = AisYf{Z), where A{s) is defined by 

^^■^^' ^^'^ " 7! U-W6 ,A^i ) (^a„(v,o + V.5:)|aa(^) j 

for s e T\Y and 

/cxp(-e-/^Z) 
•'^^^^ l^exp (-e— /3^) 



Lemma 3.5.1. Define the amplitude of the West boundary layer by (3.4-5). Let 
(si,Si+i) he a portion o/Fvi/. 

Then there exists a constant co such that the following estimates hold: 

• Far from the West end of S, i.e. for \s — sj^| > C{)5y, and for all s G 
Supp((^+(p7_i), 

\A{s)\ < C, 



(3.46) \A'is)\ < c\l + J-^^J2hyis)-y,\<s. + \cos9\Miyis)) 

(3.47) \A'\s)\ < c(^ ' ('^°^')' 



(3.48) \d^A{s)\<Ci'-'i-i^ for k £ {0,1,2}. 



• If Is- sf'l < coSy, 
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Proof. Consider a component [si_i, Si] of the West boundary. There are two 
zones on which the s derivatives of A become singular: 

• near Si-i and Si, or, more generally, near points of dQ such that \y{s) — 
yj\ — {yj 7^ 2/1 ); on such zones, the derivatives of the function ijj'^ become 
unbounded, while ■i/'i; is identically zero. 

• near s]^, the end point of E; we recall that the trace of ^^ is discontin- 
uous at this point, and that this discontinuity is lifted by the term -0''^'. 
Additionally, the s derivative of 9„V't' is also discontinuous, and this dis- 



continuity is corrected by the term ^g"''' introduced in equation (3.44 1. 
Notice however that on this zone, cos 9 is bounded away from zero. 

Away from these zones, A and its derivatives are bounded, and cos is bounded 
away from zero. Whence we now focus on the two pathological zones. Since the 
proof of the estimation on each zone is rather different, we separate the two. 

• Estimate near points such that \y{s) — j/j| ^ 1 {yj ^ yi): 
In the vicinity of such points, thanks to assumption (H3), "0? is C°° (with 
derivatives which are not uniformly bounded). We also use the following bound: 



sup 



\>^w\^)\<C- 



^^ + l 
1/4,1 + 3 



in case (H2i), 
in case (H2ii), 



so that sup^gg^j / - + -j |A^"^(s)| ^ v^/^ <^ Sy. In the rest of the proof, we write 



3.1.2 



we infer that for all s G Supp((^j ft-i): 



A for Xw- Using Lemma 

(3.49) \A{s)\ <C{1 + \Xis)\-'M{y{s))) < C. 



The estimates of the s derivatives are a little more involved. Using (3.45), we have 

A' 



A'{s)\ < C |9.0,%o(«)l + \M^'\ |5.5„0°a^(s)| 



A2 



\9ni^t\anis)\ 



Notice that on the support of ip^ (pf_^i, 



//3| 



< 



I cos 6*14/3 



< I COS ( 



thanks to the validity condition (2.8). Lemma 3.1.2 then implies, for s e Supp((/Jj fl_i), 



cost 



y' y' jei+ 



In a similar fashion, 
\A"is)\ < cf|(aM)|an(s)l + 



|(9,a„V?)|ao(s)| 



+C\\X\-'\\{did^^^)\on{s)\ + 



A2 



|(9„0°)ian(s)| 



We have 



A2 



jyl/3 



I cos 6* 1 7/3 



{\e'\' + \0"\\cose\). 
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Since cos 9 vanishes always at a higher order than 6' , we infer that on Supp(93j 'pt-i)^ 



<C, 



<C. 



I cos 61 1 7/3 

Eventually, we obtain, for s e Supp((/?^<y9^]^), 

(cos6')2 



\A"{s)\<C 



5l\\n6y\ 



• Estimates near the end point of E: 

As explained above, since the interior term is now 1/)° + '4>t., the coefficient A{s) 
changes abruptly near s]^, on an interval whose length is of order 5y. Moreover, 
since (H3) is satisfied, we have 

cos9{sY)^G, 

so that (p~(pf_i and its derivatives are 0(f) in a neighbourhood V of s^ indepen- 
dent of ly. As before, we use (3.45) and we obtain 



< 



.V3||^0 



-V^sll 



W^- 



\\Ml-^{v) 

The T/F*''°° norms of the right-hand side are to be understood in a neighbourhood 
of size 0(1) in ft of (x(sY),y{sY))- We now recall that thanks to the addition of 



^w 



the corrector term T/iff '', the amplitude A' is differentiable at s 

We now estimate W^pt + V'e||vi^'=.~ using Lemmas 
the explicit definition of ■0''^*. We infer that 

U°t+4>s\\w^^^=0{iy-^-^^) for^^e{0 
which leads to the desired estimates. 



and 



3.1.2 


3.4.1 


3.4.5 



,3}, 



together with 
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CHAPTER 4 

Proof of convergence 

The approximate solution ipapp we have buih in the previous chapter 
(i) satisfies exactly the boundary conditions 

Wapp\QQ^ — \^nWapp)\d^ ^ t) 

by definition of the boundary layer terms; 
(ii) has H^ regularity, but with non uniform bounds on the derivatives (see for 

instance the definition of the interior term near East corners and surface 

discontinuities); 
(iii) but does not satisfy the Munk equation, only the approximate version 

dxtpapp - V^^i^app =T + 5t, 

where 6t comes from both the interior and boundary terms. 

The goal of this final chapter is to estimate the different contributions to this 
error term 6t, and to check that it is admissible in the sense of Definition |1.3.2| 
At this stage, there are no more conceptual difficulties: all the proofs are based on 
technical computations. 



4.1. Remainders stemming from the interior term t/)™* = t/jJ" + ?/)'* 

Going back to the construction of the interior term, we see that its contribution 
^Tint to the remainder is made of two kinds of terms. 

• Because of the truncation, the solution ip^ to the transport equation sat- 
isfies the Munk equation with source terms r(l — Xi^) and —I'A'^ip^, which 
are admissible remainders according to Lemma |3.1.1[ 

• The lifting term which restores the H^ regularity near discontinuity sur- 
faces also introduces correctors. Part of them are considered as source 
terms for the singular surface boundary layers and will be dealt with in 
the next section. The other ones will be proved to be admissible remain- 
ders as stated in Lemma [2.4.31 

4.1.1. Error terms due to the truncation Xi^- 

Let us first prove Lemma [3.1.1[ i.e. the following estimates : 
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• We start by localizing the problem around {xi, yi) thanks to a partition of unity: 

for i G /+, let f, G Cg°(R2) such that 

^ /,; = 1 on VI, 

and 

• for all i G /+, /i = 1 in a neighbourhood of {xi,yi), 

• for all i ^ j, d{{xi,yi), Supp /j) > 0. 

Then, according to the definition of Xf j for ^ small enough, 






Thus, defining 'ipf by 

SxV'j' = fiTXz {x-Xi,y-yi), 

we have ipt — '^''Pi ^ ^-i^d it suffices to prove that under the assumptions of Lemma 

[MS 

\\fir[Xtix~x^,y-yi)- 1]||h-2(o) =o{i^^'^), 

|kAV°||//-(a±) = H|AV."||L^a^±) - o(i.5/8). 

In order not to burden the notation, we drop the indices i and we shift the 
origin of the arc-length parametrization so that Si = 0. We also shift the origin of 
the axes so that (z^, yi) = (0, 0), and we rename respectively t and x^ the functions 
fiT and Xi(a; — Xi, y — j/i) (which have the same regularity as the original functions). 

Note that, by definition of x^, we have the obvious estimates 

K'd;Hi - x„)iu. < Cm« ((^)" [^)" (Sis,)"^) ■ 

• We now give an upper-bound for the source term t(1 — Xi^) in terms of J^, dy. 
Without loss of generality, we can assume that s = is a North boundary point of 
r^. Define, for (a;,y) G fi, 

/y ry' 

dy' / dz{xu - l){z). 
-OO J —OQ 

so that 

dy(f>iy = Xi^ - ^, (f>i^{x,y) = ioT {x,y) en, d{{x,y),dQ) ::^ Sx,5y. 
Since the support of 1 — Xi^ is included in the rectangle 

[^S-{l + \\nSy\-'),St(l + \\n6y\-^)] x [-Sy,6y], 

we have 

(l>^{x, y)^0iix>6+{l + \ \ndy\-'^) or x < -6-{l + \ In^^Ti), 

and (puix, J/) = if y < -Sy. 
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Figure 4.1. cos 9 has a local minimum at s = 

Notice also that 

\Mx,y)\<C{6l + \y\Sy) V(x,2;) e a 
Therefore, parametrizing locally dfl by a graph {x,yN{x)), we obtain 

\\M\hin)<C f ^ r \sl + Sy) dy dx. 

In order to bound the right-hand side of the above inequality, we have to distinguish 
between two cases, depending on whether cos 9 has a local extremum at s = 0. 

• No local extremum at s = 0: either cos 9 changes sign across s = 0, or 
cos 9 is identically zero in a neighbourhood on the left or on the right of 
s = 0. Notice that in this case, xsiy) is only defined on a neighbourhood 
of the left of y = 0, and thus we have 5~ = 6^ =: Sx- Moreover the 
function yN{x) remains non-positive for |a::| < 5x{^ + |ln(5y|^^). In this 
case, we simply obtain 



mw^n) < cs%. 



• Local extremum at s = (see Figure 4.1) : we have 



M\hm < CSX + C f iSt\yNix)\ + 6l\yNix)\^) dx. 

We now estimate the right-hand side using the formulas in Appendix A : 
yN{x) ^ Cx""*"^ in case (i), 

j/Affa;) ^ Cx^ exp { ~- — r ) in case (ii). 
V \x\ I 



As a consequence, 

(,5-)"+2 in case (i), 

in case (ii) 



/■o r (5-)»+2 

/ \yN{x)\dx<C\ (X-\i^,,^( g A 
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/ 



and 
We infer that in all cases, 



3a 



+ |ln5J-i) 



in case (i), 
in case (ii) 



\\Mhin)<CSl{5++S-). 
Therefore, by definition of 6^ and 6y, we get 

since in the worse case 

Sx = max((5^ , S~ ) 



a 4a 



• It remains then to estimate the viscous term i/A^t/j" in H '^ (il^), that is to 
estimate vAip^ in L^(51^). We indeed recall that, by Proposition 2.4.2 

We then start from the identity 



A^o = 5.(rx.)-52 



'E{y) 



(4.3) 

(4.4) 
(4.5) 



XE{y) 



= XudxT + rdxXu - / (Xi^dyT + 2dyTdyXv + rdyX^) 

J X 

~'^x'j^{y){xudyT + TdyXi,)\x=XEiy) ~ XE{y){'^Xu){xE{v),y) 

-WEiy)? iXudxT + TdxXu)\x=XE(v) ■ 



Using (4.1 1, we can easily check that 



llQlU^ < max (l + Sl/'iS^)-'^'\lnSy\ + S^S^^^^) . 

Note that we use here in a crucial way the fact that there is no discontinuity line 
in the domains fl^ : the singularity is localized in a rectangle of width Sx- 



As for (4.4), ( |4.5| ), we use (3.11) together with the formulas of Appendix A: 
• Case (i): we deduce that 

\ 1/2 

WEiy)\^dy] +C 



1(4.4)1 



C 
< 



\y\>Sy/2 



1/2 



WEiy)\'dy 



\y\>^yf2 



< CSy 



C6y 



1(4.5)1 



C 
< max -rr: 
- ± 5± 



< CSy 



1/2 



\y\>Sy/-2 



WE{y)\'dy 



CSy ^<"+^' . 



Hence the choice Sy = z^^/^l ln(z^)|^/^(ln | ln(i^)|)"^, with ^ > arbitrary, 
ensures that St is an admissible remainder. Notice that any choice of the 
type Sy = ly'^, with 7 e ( ^5^+6) ' 4(3n+ 1 ) ) ' ^^^° works in this case. 
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Case (ii): similarly, we have 



» ^ id 



1/2 



y|>c5„ 2/^(ln|2/|) 



4C^y 



< '5;3/'|ln|'5,||- 



C 



1 



^' /M>c.,ynin|y|)4'^^ 



1/2 



1/2 



± 5.j; \J\y\>c5y 
1-3 



(4.5 II2 < max^:^ / dy 



Thus, with the choice Sy = j/i/"*! ln(z^)|i/^(ln | ln(i^)|)-'5, we infer that St 
is an admissible remainder. 



This completes the proof of Lemma |3.1.1| 

4.1.2. Error terms due to the lifting term ?/;''^*. We give the rest of 
the proof of Lemma |2.4.3[ Since we have already checked that all the terms are 
well-defined, there only remains to estimate the size of the remainders coming from 
ip^^^^. More precisely, we prove that 



(4.6) 
where 



'lift 



4u 



llriJU.(^) = 0(^1/8)^ WrUn-^n) = 0(1^'/'), 



'"x"^ (^) h"{x) + 2.'/\+x" (^) («"(x) + b"{x){y - y,)) 



rL = rix. - 1) + ^X+X (^) («*'H^) + b^'Hx){y - Vi)) 
We rely on Lemma |B.1| in Appendix B, which yields 

W"^' ("9?) ("" + ^"(y - 2;i)) + 4.^/V (^ 






b" 



L2(o) 



As for Tj^jfj, we have already proved (see Lemma 3.1.1 ) that t{xi, — 1) = o(j^^/^) in 
H^^{n). And we have clearly 



^x+x (jjTjr) («^'n^) + b^'H^){y vi)) 

{a"{x)+b"{x){y-y^)) 



^l ( '^X+X ( ^^^^ 



so that 



\\rU\H-^in) - 0(^5/8) + o (^9/«(|la"|U. + ^i/^||&"|U.)) = o{v''^). 



Gathering all the terms, we obtain (4.6) 
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4.2. Remainders coming from the boundary terms 



Following the construction of the approximate solution, we now turn to the 
remainders coming from boundary layer terms. This step is more technical since 
boundary terms are defined through rescaled curvilinear coordinates, so that esti- 
mating the bilaplacian requires tedious computations. 

Actually most of the terms will be estimated in H^^ norm, so that we will need 
to compute only one iteration of the laplacian. 

4.2.1. Laplacian in curvilinear coordinates. In order to get the boundary 
layer equation (2.3), we have kept only the leading order term of the bilaplacian. 



i.e. the fourth derivative with respect to the scaled normal variable Z. Rewriting 
this boundary layer equation as a Munk equation with remainder terms, we have 
to estimate the lower order terms of the bilaplacian, which can be done thanks to 
the following 

Proposition 4.2.1. Denote as previously by {z,s) the local coordinates on a 
tubular neighbourhood of the boundary. Let f = F{s, A(s)z) be a smooth function 
with Supp / C Supp xq ■ 

Then, we have the following error estimate 

|A/ - d,J\ < C {\dssF\ + {y/Xf\dzzF\ + \y/X\ \Zd,zF\) 



C 



\Zd..F\ 



\'e" 



\Z^dy,F\ 



\'' 



\ZdzF\ + \\e'\\dzF\ 



In particular, 

' if X does not depend on s 



(4.7) 



|A/ - d,J\ < C \ds,F 



A 



\zdsF\ + \xe'\\dzF\ 



• if F is an exponential profile with respect to Z , namely if F{s, Z) — 
v4(s) exp(— aZ), with a G C a constant of order one with positive real 
part. 



(4.8) |A/-a,,/| < C[\dssA\+[ - + - ]\dsA\ 



+C 





( 


0" 




\' 


,sA\ + 1 


T 


+ 


J 


\'0" 




A" 






+ 




+ lAe*' 


X^ 




A 







-cZ 



\A\e-'^ 



where c — ^{a)/2. 



Proof. We start from the expression of the Laplacian in curvilinear coordi- 
nates 



A = 



1 



d^{{l + z9')d^) + d4^^d. 



l + z6' 

Given the specific form of the profile / = F{s, A(s)z), we thus have 

Xe' „ „ z6" 



A/ = X'dzzF 



(4.9) 



l + zB' 



dzF- 



{i + ze'f 



-^(dsF + X'zdzF) 



1 



[l + zO') 



— [dssF + 2X'zdszF + X"zdzF + [X'zfdzzF] 
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Then, using the faet that on the support of xo the jacobian of the change of variables 
is uniformly bounded from above and from below 

0<^<l + ze' <C, 

we immediately get the first, general, estimate. 

For North, South and discontinuity layers, we always choose 

A = Z/-1/4 

so that A' = A" = 0. The remainder consists then only on three terms. 

For East and West boundary layers, / is a combination of exponential profiles. 
Therefore, multiplying by Z or differentiating with respect to Z does not change 
the size of the error term. 

D 

In order to get the expected estimates on the remainders, we first note that 
A^ = {X^dzz)^ + A(A ~ X^dzz) + (A - \^dzz){\^dzz) 

- {>?dzzf + (A + \^dzz){^ - >?dzz) + [A - \^dzz,>?dzz] 

• The first term of the right-hand side has been considered in the equation 
for the boundary layers, so it does not introduce any remainder. 

• For the second term, we expect that it can be estimated in H~^ using 
the previous proposition together with the controls on the boundary layer 
profiles and their derivatives up to second order. 

• The difficulty comes then from the third term, which is a commutator 
between two second-order differential operators. Roughly speaking, it 
should be a third-order differential operator, so that we could estimate 
the corresponding remainder in L^ in terms of the third-order derivative 
of the boundary layer profiles. 

However we will not proceed exactly in that way, first of all because the com- 
mutator [A — \^ d^ z T '^ (^^ z\ involves too many terms so that the computations 
would be very tedious, and overall because some terms are not small in L^ and 
need to be estimated in H~^ . The proof will actually be adapted to the case to 
be considered: recall indeed that for North, South and discontinuity layers, A does 
not depend on s, whereas for East and West boundary layers, F is an exponential 
profile with respect to Z. 

4.2.2. Error terms associated with North and South layers. 

In order to simplify the presentation, we focus on a South boundary layer, as 
we did in the preceding chapter. Of course, the case of North boundary layers is 
strictly identical. We denote by Si_i, si the end points of the South boundary under 



consideration, and we introduce a truncation function 7^-1, i as in (3.36). 



Lemma 4.2.2. Lettps be the solution to (3.271. 
Then 

with 

||ri|U.(o)=o(^i/«), \\rnH-Hn)=o{,.'/'). 
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Proof of Lemma [4.2.21 Since ips satisfies (3.271, there are several kinds of 
terms in 



• All terms in which at least one derivative of Xo or 7i-i.i occurs are small: 
indeed, we have one derivative with respect to z less acting on ips, so that we gain 
a power i/^''*. Typically 



I'dti'^sXo) - '^Xodt-ips 



■Y,Ciiy^'d'z^s{-s,Xz)dt ^xo 



^4-J, 



j=0 



(4.10) = Oii^^'/'UshiHl) ^r, L\n) 

+0{,yy''\\dzyJs\\LlJ'^^H-\n). 
In the same way, 

•05(5, '^~^^'^z)dxixo{z)%-i.i{s)) 



^s{-^,y-^'^z) 



smt/ls 



^ Qids^i-i^r{s)xs){z) +cos6l(s)7,,_i,i(s)92Xo(z) 
Using the estimates of paragraph |3.3.3[ we infer that 



ipsis,!^ ^^^z) cos6{s)ji-i,i{s)dzXo{z) 



L2(o) 



0{v^/^\lniy\-^). 



Moreover, on the support oidgji^i.i, either V's is zero or H^sHloo^^.^ ) = 0{\ Inv ^). 
We deduce that 

i;s{s,>^-'^''z)smd{s)ds^^^i4s)xa{z) = 0(z/i/«| Ini^ri), 

L^ (O) 

SO that eventually 

(4.11) ||^5(s, iy-'^'z)d,{xoizH-idsmLHn) = 0(i^i/«| Inz.]-!). 

• The terms stemming from the d^ derivative applied to ^5 can be split into a 
remainder term coming from the jacobian 



zd smt' f ^fi I f -1/4 N 

1 + z0' ^o(zni-iAs)OsVs{S: V ' z), 



and 



XQ{z)l^-lAs) (y-^'^ cosedzi^s{s,v-^'^z) + smedsiJs{s,y'^'^z)) , 



which will simplify with the term v\ XG{z)'yi-i,i{s){d\)'4!s{s^ v ^'*z) coming from 
the bilaplacian. 



(4.12) 



From the moment estimate in (|3.3l|), we thus deduce that 



ze'sinO^ , .p. , , „i/4 N 

Y^p^7i-i,»(5)Xo(2)OsV'S'(s,i^ ' z) 



L2 



• The most technical part is the control of the remainder terms coming from the 
bilaplacian. For the sake of simplicity, instead of commuting the whole laplacian 
with A^9|, we will first get rid of the corrections coming from the jacobian, then 
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commute the remaining part of the laplacian. More precisely, we start from (4.9) 
with a constant A 

= A(A - \^d\) + (A - \^d\)\^d\ 



A(A - X^dzz) + 



\B' 



-.d. 



ra. 



l + zd'^^ {l + zO'f " {l + z6')'^ 
A(A - X^dzz) + X^dzzixe'dz - zB"ds + d,\ - 2Xd"dz, 



dss 1 A^a. 



zz 



xzio'f. _ £^014)!^, ^ (l±^^tla. ] x^d 



l + zO' 



{l + zO'f 



[l + zO'f 



zz- 



We have, using the regularity and moment estimates (3.311 
i^(A - X^dzz)^s 

xe' x-^ze" ^ V ^ 



1 + ^6" 

We now use the following commutator estimates, of which we postpone the proof: 

Lemma 4.2.3. Let p e H'^{Q.), with support in {d{x,dQ) < S}, and let f e 
L^{dnxR+). We setX:^v-^l^. Then 

2 



\p{s,z){X^dl)f{s,Xz)\\^_ 



'■m 



< 



C\X\-'/'\\f\\LHon.u^)Y.\\9'p\\^- 



in), 



fc=0 



\\p{s, 2)A(/(s, Az))||^_.(^) < C\X\-'^^f\\man.B.^)\\p\\w^.^in)- 
The second commutator estimate of Lemma |4.2.3| implies 

(4.13) \\xo{zn^iA^>M^ - )^^dzz))Ms, Az)IU-.(o) = o(^'/'). 
In the same way, we have 

vXO'dztps - i^z0"dsijjs + i^dssi^s 
= 0(^.^.-1/^1/8)^2 + 0(1^^.1/4^^1^.1/8)^2 + 0(z.i.-i/2| In i.|-ii.i/8)i2 
so that, using the first estimate of Lemma [4. 2. 3( 

(4.14) i^Xo{zn.-iAs)^^dzz{XO'dz - zO'ds + dss)tps = o(^'/')h-= ■ 

For the corrections coming from the jacobian, we use the first estimate of the 
commutator lemma 14.2.31 

fXzie')\ Z0"((l + Z0')'-l))o (l + ^^')'-lo \ ^2^ , 



H-2(n) 



= o(:.i/8)||zazVsllL^ 

+ 0(;.l/8t.l/2)||Z25,^5||L^ 

+ o(^l/8^l/4)||z52^5L^ 



'\2 



1 + ZI 



^Xo7j-i,i 



^^■^(a) 



'((1 + Z0') -1) 



z(l + z6'')3 

e"i{i + ze'f -1) 
z{i + ze'y 



■XOli-lA 



-Xoli~i,i 



w^.^(n) 



W2,oo(Q) 
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which imphes that 
(4.15) 

The only remaining term is 
(4.16) 2v\H"xo{z)l^-iAs))dzs^s = 0{vv-^'H-^v^'^)l2. 



Combining estimates (4.10)- (4. 16 1 shows that the remainders coming from the 



South and North boundary layers are admissible in the sense of Definition 1.3.2 D 
We now prove the commutator estimates: 

Proof of Lemma [4.2.31 

• Let us first recall that \dz ■ iI*(R+) — > iJ''^^(R+) is a bounded operator 
(uniformly wih respect to v) for all s. 

Let ifi G H^{fl) be an arbitrary test function such that (p\gfi — dnf\dn = 0. 
Then, using the Cauchy-Schwarz inequality, 

p(s, z){\^dl)f{s, \z)Lp{s, z){l + z6') ds dz 
n 

/(s, Az)af (/9(s, z)ip{s, z)(l + z9')) ds dz 
n 

< C|A|-i/2||/|U.(aoxR+) E \\d'AP^)\\L^n) 

2 
k=0 

By definition of iJ-2(^)^ 

||p(s, z)(A2(9|)/(s, Az)||^_,(j^^ 
= sup / p(s, z)(A^9|)/(s, Az)(^(s, 2:)(1 + 26*') ds dz, 

llvllH2(n,<l 

from which we deduce the first commutator estimate. 

• Concerning the second inequality, we merely develop the commutator [A,p]. 
We have 

p(s,z)A/(,s,Az)-A(p(,s,2)/(s,Az)) = -2VpV(/(s, Az)) - Ap/(s, Az) 

= -2div(/(s,Az)Vp) + Ap/(s,Az). 

The estimate follows. 

D 
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4.2.3. Error terms associated with East and West boundary layers. 

We recall that the boundary layer term on the East or West boundary [s^, s^+i] 
takes the form 

^E,w{s,Z) = {^+ip-^^){s)A{sYf{Z), 
where 

• A{s) — — A^^(s)9„'i/;"|p (s) and f{Z) = exp(— Z) on East coasts, 

• yl is the matrix defined by (3.45) and 

exp (-e"/3Z) 
^exp (-e-*"/3z) 

on West coasts. 



liZ) 



with 



Lemma 4.2.4. Let A^f he defined by the expressions above and (3.45) 
Then 



Proof of Lemma [4.2.41 Throughout the proof, we write A instead of \e,w^ 
and we set Lpi = ip'l if'^,-^^. Since / is an exponential profile, we will use the second 



part of Proposition 4.2.1 to estimate the remainder terms in the bilaplacian. In a 
first step, we bound the error terms r\^ yy and r|; ^y by expressions involving v, 
and j4, and in a second step, we prove that these bounds yield the desired estimates. 
First step: estimates in terms of the amplitude A 

We claim that for all i such that (s^, s^+i) C Ty/ n F^, 

a - u/S.^){^,{s)A{s)f{\z)) = t\ + r^, 
where 



+C,y^^''\\A\\L^^an), 



(4.17) 
and 

(4.18) +C ( |{^s{v^A)f 



1/3 \ 1/2 

ds\ 

5/3 \ 1/2 

ds\ 



+C( A(s)2(0')'t^'/'|cos^r"/3ds 

\JSupp ipi 

+C|| A||o,i/5/6 + Ci^-^ exp(-C/z/i/4). 
We recall the expression of dx in curvilinear coordinates: 




1/2 



d_ 
dx 



d sin 6 d 



dz l + zO'ds 
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Throughout the proof, wc use the fohowing estimate: 
(4.19) yseSuppipi, Ab,vi/(s) > Czy-i/4. 



This is a consequence of assumption (H2) together with Lemma 3.2.1 
We now estimate 

[d^ - z/A^] {ip,{s)A{s) exp(-A(s)z)xo(^)) • 

• We start by commuting the differential operator dx — vt^ with the multiph- 
cation by xo- Since Xo(-z) = 1 for z in a neighbourhood of zero, all terms in which 
at least one derivative of xo appears are exponentially small: indeed, for all fc > 1, 
and for aG {1, e^^^'^/^i.^ 

|i9^Xo(-2) exp (— q;A(s)z)| < Cfcexp(— Ct^""'"' ) V(s,z) G Supp(^i x Suppxo- 
We also have (see Lemmas 3.1. 2[ 3.4.1 3.4.5) 

|a^a,^(¥>i(s)A(s)exp(-aA(s)z))| < Ci^"-^ exp(-cA(s)z)), 
for some positive constants c^C^K . Hence we infer that 

[dx - v/\^] (^i(s)yl(s)exp(-A(s)z)xo(2:)) 
= Xo{z)[dx - v/^^] {^^{s)A{s)c^p{-\{s)z)) 
+0{v-^ ciip{-Cv-^/^)) in L^{n). 

• The terms stemming from 9j.((^i(s)A(s) exp(— A(s)z)) can be split into 

\coseip,{s)A{s){dzf){\z), 

which will simplify with the term v\'^Lpi{s)A{s){dy)f{\z) coming from the bilapla- 
cian, and a remainder term 

sin " 



(4.20) 



1 



zt 



- [ds{f^A)f{\z) + A{s)^,{s)zX'{dz!){\z)] 



In order to replace the jacobian term by 1 in the expression above, we use the 
following estimates, of which we postpone the proof: 

Lemma 4.2.5. Let I C Ve^^w be a closed interval, and let A G L°°{I, C) such 
that inf/ 5R(A) > 0. 

Let 4> G L°°(il) such that Supp C /. Then 

1/2 

ms)xoiz) exp(-Az)||^.(^) < C ^ 0(s)^(5RA(s))- ' 
and 



< C 



s)xo(2)exp(-Az)|l^_2(f2) 

1/2 

is)''i^Xis))-^ds' 



C||(/.||ooCxp(-Cinf3fiA). 



Therefore, replacing {l + z9') -'^ by 1 in (4.20) yields an error term whose square 



L norm is bounded by 



C I 9'^^ 

/Suppipi |C0S( 



\ds{Acp,)\' 



ql2 



COS^ 



-}A{s)?] ds 



< 



C^\ln,.\\m^,A)\\l^^g,,^+C,.'^/^\\A\\l^^g^y 
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Indeed, by definition of Lpi, 



Supp(pj C (s+,s^+i) 



and using assumption (H2), it can be easily proved that for all i such that (s^, Si+i) 
is a western or eastern boundary, 



e'^\co&e\-^ ^0{\\nv\), 



'icOS^r^^ 0^2/7)^ 



Then, using Lemma [4. 2. 5[ we infer that 

||sin6'xo(z)as(¥'iA)/(Az)||^_2(n) 

5/3 



< c / {ds{v.,A)Y 

Jon 



cos{e{s)) 



ds + 0(z/-^exp(-Cz^-i/4)) 



and 



\sm0xo{z)A{s)Lpi{s)zX'f{Xz)\\^_2(^^) 



< CI A(s)20'2z.5/3|cos0|-"/3rf<, + o(j,-ifexp(-Ci^-i/4)). 

' Supp ^i 



We now develop the bilaplacian of {ipiA){s)f{Xz), using Proposition 4.2.1 



We bound all the derivatives of by a positive constant. Moreover, thanks to the 
compatibility condition (2.8) and the definition of iff, it can be checked that for 
s g Supp ifi, 



< 1, 



A" 
A2 



«1, 



(A 



/\2 



A3 



<1. 



Therefore we obtain 



UA-X^^zz){v^A)is)f{Xz)\ 



< c \{Aip,r\ + 



|(A(^,)'l + |A^<^»| exp(-c|A|z) 



We now evaluate the right-hand side of the above inequality in L^(fi) (after multi- 
plication by xq)- By Lemma 4.2. 5[ we have 

\\{XA^,){s)exp(-c\X\z)xo{z)\\mn) < C\\A^A'^^\\l^ < C\\A\\^iy-'^^ 
Therefore 

A{iAip,)is)fiXz)) = iAip,)is)Xy"iXz) 

+ O {\\A\\l^i.-'^' + ||(A^,)'A'A-3/'|U2 + ||(A^,)"A-'/'||l^)^^ 
Commuting once again xo and A, we infer that 
,^Xo{z)A^ i{A^,){s)f{Xz)) = :yxoiz)A ((A^,)(s)AV"(Az)) + 0{expi~C,.'/*))L2 
+ 0(P||oot^'/' + t^||(A^.)'A'A-3/2||^, + iy\\{A^,)"X-'/^h2)H-2 
Notice that on Supp(y9,i, 

I.|A'|A-3/2 « A-5/2, 
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SO that 



y\\{A^,)'>^X-^'^\ 




{^s{v^A)f 



cos{9{s)) 



5/3 \ 1/2^ 

ds\ 



It remains then to deal with the first term i/xoi^)^ ((^<<5i)(s)A2/"(Az)). We 
develop the laplacian one more time. Since / is an exponential, we can use the 
computations above and simply replace A by AX^. Remember also that thanks to 



(2.8), on Supptpi, we have 



Therefore 



A''a-i/2 « x'/\ 
|A"|Ai/2 « A5/2. 



z.Xo(^)A((A^,)(s)AV"(Az)) 
= i.xo{z){A^:){s)X'f^'HXz) + j^f^^iA^.rx'nXz) 

+,yO{\\Aip.,X'/^L2 + \\iAip.,yX'X'^^\\L2) in L^{n). 

We estimate the second term of the right-hand side in H~^{Vt). This is the only 
time where we really need to commute the jacobian terms with the z derivatives, 
as explained at the end of paragraph |4.2.1[ Using Lemma 4.2.5[ we have 

Xq{z) 



{i + z9'y 



r(A^O"A'/"(A^) 
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Xo{z) 



(Aipi)"f{Xz) 



dz^ \{l + z9')^ 

,. Xo(z)9' ,,,,, ... y, „ Xo(z)9' , .,. w, 



L2(Q) 



+0(||(A^,)"A'/'||L^exp(-C/^.i/4))^ 



Notice also that on Supp tpi 



i.A'Ai/2 



<c\e'\ 



cost 



1/2 



SO that 



\{A^,yX'X'/'\\L2 < C||(v4^,)'lloo^'/'|ln:.|i/^ 



Gathering all the terms, we obtain (4.17), (4.18). 

Second step: quantitative bounds for rj, rf: 

We now use the definition of A together with Lemma [XTT2] and Lemma [XO] in 
order to estimate the right-hand sides of (4.17), (4.18). Notice that the estimates 
on West coasts are always more singular than the ones on East coasts: indeed, (pi 
has bounded derivatives on East coasts, and unbounded on West coasts. Moreover, 
the estimates of Lemma [3]T]2] show that A is always larger on West coasts, together 
with its derivatives. Therefore, we focus on a West part of the boundary in the 
following, i.e. we assume that (si,Si+i) C F^. 
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Following Lcmina 3.5.1 we treat separately the estimates in the vicinity of Si 
and in the vicinity of points s such that cos 9 = or y{s) = yj for some j G /+. In 
the rest of the domain, the estimates on rj and rf are easily proved. 

• Estimate near points such that \y{s) — j/j| <C 1 {yj 7^ yi): 
Let us first recall (see Lemma 3.3.21 

||a,V»IU<C7j/-'=/^ fee {0,1, 2}. 

Therefore 

Wrlh^ < C [i^'^'WAlU + i^'/'l ln^.|i/2p'|U + i^'^'\\A"\\ 

Inequalities ( |3.49[ ), ( |3.46[ ), ( |3.47[ ) yield respectively 

(4.21) ||A||oo-0(l), \\A'\\^ ^ {S-'\\ndy\-'), \\A"\\^ ^ 0{d;'\ln5yr'), 

so that llj'j^llLa = o(i/^/^). We now address the bound of rf. Using once again 
( |4.21[ ), we have 

1/3 \ 1/2 
ds\ 



\Jdn 



cos{e{s)) 



< Cv'l'' (||^"|U + ll^'llooP'lIco + P"||oo) U ISupp^J COS^rl/^") 

< Ci/2/3|lni.|-ir/"lsupp^Jcos0|-i 



1/2 



Using the formulas in Appendix A together with the definition of 5^,5,-, j^, it can 
be proved (treating cases (i) and (ii) in assumption (112) separately) that 



(4.22) 

We infer that 



/ |cos6l(s)r''ds = 0(z/-''/'^|lni/r(2+37/2))^ 

-/ Supp ipi 



\Jdi1 



cos{e{s)) 



1/3 \ 1/2 

ds = o{iy^/^). 



In a similar way, we have 



a'^2| 



Supp ipi 



,S0rll/3^O(^^-5/12|ij^^|-l 



/2 



SO that 



1/2 



A(s)2(0')'i^'/'l cos0r"/3 ds = 0(iyt-&| lni.|-i/4) ^ ^(iyS/S). 

'Suppi/Pi / 

We now tackle the term 

5/3 
A'is)id,ip,{s)f 



(4.23) 



on 



cos(6'(s)) 



ds. 



By definition, the support of dgipi is located in neighbourhoods of Si and Si+i, so 
that we can use assumption (H2). We distinguish between the cases (i) and (ii). 
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At this point, assumption (3.24 1 is necessary. In a neighbourhood of s = Si, in case 
(i), we have 

\dsipt{s)\<Ci^-^l _^ _^, 



SO that (4.23) is bounded by 

5 1 i 1 l"i &n\ / 1071 + 11 \ , -/., 



In case (ii), using Lemmas 3.2.1 and 3.3.2 we infer that 

In I In z^ I 



Is±-a±| 



C- 



while 



(lni.)2 ' 
|cos0r^/3<C:/-5/i2|lnz.|-9/2. 



We deduce that the contribution of (4.23) is bounded by 

|sf -a±|-"z/^/3 / I cos 01-5/3^0(1/5/^1 In !.ri/2(ln I In i.|)-2)=o(z/^/4). 
There only remains to check that 



./Supp t^i 



cos{e{s)) 



5/3 



ds = 0(l/5/4). 



We use inequahty (3.46). We have 

f (^'(s))2|cOS0|-5/3d 

cf \cos9\-^/^ds + C f 



cos 



llppl^i 

c 



Supp '^i 



0\"^M{y{s)) 



\y(s)-yj\<&y 



Using (4.22), it can be checked that the first term is o(i/~5/i2)^ There remains to 
estimate the second and third term of the right-hand side. We have to distinguish 
between two cases: 

• If I cos 6*1 ^ 1 in a neighbourhood of the point where \y{s) — yj\ ^ 1, then 
we can use assumption (H2) and the estimates of Appendix A. We infer 
that 

/ Icos^l^/^ ^ l\yis)^y,\<8, < C6l^'\\ndy\-'^\ 



cos 



0|l/3X(zj(s)) < CC'^3|1j^j^|-10/3 



/Supp ipi 

• If COS 6 is bounded away from zero at the point where y{s) = yj for some 
s G Supp (fii, j e /+, then we can use the change of variables s — > y in the 
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integrals. Since dy{s)/ds = —cos9{s), the jacobian is not singular, and 
we have 

f \cose\'/^Y.hvis)-yA<Syds<C f ldy<CSy 

COS0|l/3ls 



' + 



inf.e/, \y{s) - y.min \y{s) - y,\r H l^»<2|.(«)-..l<i 



< c \y\-^My\)-' 

JSy/2<\y\<l 

In both cases, we infer that 

1/5/3 f (A'(s))2|cos0|^5/3^s = o(iy5/4). 

-' Supp ipi 

This concludes the estimates in the neighbourhood of a point s e dft such that 
COS0 vanishes or y{s) ~ yj for some j € I+. 

• Estimates near the end p oint of E: 

Using once again Lemma 3.5.1 we infer that ||rj^||L2 is bounded by 

while ||rf 11^-2 is bounded by 

This concludes the proof. D 

There remains to prove the Lemma giving the sizes of East and West boundary 
layer terms in L^ and H^^: 

Proof of Lemma 14.2.51 We begin with the L^ estimate. We recall that the 
jacobian of the change of variables {x,y) — >■ {s,z) is equal to (1 + z9')~^, and is 
thus bounded in L°° . As a consequence, 

\\(t){s)xo{z)exp{~Xz)\\l2(^n) 
-s 

i){sfxo{zf exp(~25R(A)z)(l + zO'y^dsdz 
dnJo 

2 



< C / 0M / exp(-23ff(A)z) J dsdz 

< C f (t)^i^Xr\ 
As for the H^^ estimate, we have 

^2 

<?^(s)xo(-2)exp(-Az) = — (0(s)xo(2)(A(s)) 2exp(-Az)) 

+20(s)x^(z)(A)-iexp(-A^) 
-Hs)Xoiz)exp{~Xz). 
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The last two terms are supported in {S/2 < z < 6}, hence they are exponentially 
small in L^. More precisely, they are bounded in L°°(ri) by 

C||</.||ooexp(-(5inf3fi(A)/2). 
As for the first term, if C G Hg{n) is an arbitrary test function, then 
0^ (0('S)xo(^)(A(s)) ^ exp(-Az)) (^{x, y) dx dy 

('/'(s)Xo(z)(A(s)) "^ cyip{-\z)) C,{s, z)——^^ds dz 



1 1 Jo 9z^ 



d^ 



uo 



(s)XoW(A)-^Gxp(-Az)— ^CGs,z)-^^^g 



ds dz 



1/2 



< ^ll^ll^^(f^) lyj^ 0(s)^|Ar^exp(-2$R(A)z)dsdz 



< c\\aHHn)(^l^cb'mr'^ 



1/2 



By definition of the H norm, we infer the estimate of Lemma 4.2.5 



n 

4.2.4. Error terms associated with discontinuity layers. 

Since we have already proved that the error terms associated with ip^^"^ are 
admissible in the sense of Definition 1.3.2 (see Lemma 3.4.5 1, there only remains to 
prove that 



/Q A2^ I ;2J I y~yi\ ( y-yi 
(a, - i/A ) u [x, -—j^ ] X 



(5r"" X, 



1/1/4 J ^\ Sy J J V ' !^l/4 



y-yi 



'E: 



where 

Note that the computations here are much simpler since all the formulas are given 
in cartesian coordinates. 

By definition of ip^ , we have 



y-yi 



y-yi\ .^lift / ^ y-y i 
y-yi 



(4.24) = (x-1) 
(4.25) 

(4.26) 



y-yi 



^-"" -'^ 



y-yi 

1/1/4 



MH + d'M ( V^'^ ( X: 



A/i 



y-yi 



^v 



3=0 



- < x,'-^) /'-^^ ('-^ 



The term (4.24 1 is supported in 

{\y~yi\<'^'/'}n{\y-y^\>Sy}. 



3.4.1 



Since Sy ^ u ' , (4.24) is zero for i/ small enough. Moreover, using Lemma 

we infer that (4.25) is 0(i/i/^^/^^j/^/*) in H^^{fl). Eventually, using the same 
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type of commutations as in the proof of Lemma 4.2.2 we prove that (4.26) is 
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Appendix A: Equivalents for the coordinates of boundary points near 

horizontal parts 

In the vicinity of a cancellation point sq G OTe, the East part of the boundary 
can be described as the graph of some function y i— >■ XE{y)- Assumption (H2) 
provides then asymptotic expansions of the function xe ■ 

Let (a;(s),y(s)) be the coordinates of the point with arc-length s on 551. By 
definition, if cos(0(s)) > 0, 

x{s) ^XE{y{s)). 
Moreover, 

d (x{s)\ _ ( sin(6i(s)) 



^^■^^ ds\y{s))-\-cos{e{s)) 

and 

x'e{v{s)) = -tan(6l(s)), 



COS3(0(S))' 

We set the origin of the axes so that (a;(so),2/(so)) — (0,0). 
In case (H2i), we have 

y{s)^C{s-Sor+'^Cx{sr+\ 

and thus 

XE{v)-^C\y\^^, 

(A.3) x'E{y)-C'\y\^-, 

x'^{y) ^ C"\y\^^'\ 

In case (H2ii), using the same kind of calculations as above, we infer that as 
y vanishes 

C 



XE{y) 



(A.4) x'E(y) 



ln|2/| 
C 



j/(ln|y|)2' 

y\\n\y\y 
We have used the following fact: 

( " ^ J ' {s- sof / a , 

exp as ~ exp as s — > sq ■ 

so \ s' - soj a V s - So . 
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Appendix B : Estimates on the coefficients a and b. 

Lemma B.l. Let 

Xmin := inf{a; G R, 3y G R, {x, y) G Vt}, 
Xmax ■■= supjx G R, 3j/ G R, {x,y) G Q.} 
Then for aZ/ fc G {1, • • • , 4}, 

I!«"IIl^(.„.„,.™„.) + ^'/'II&"IIl^(.„„„,.™„.) - o(^'/')- 



Proof. We recall that a and & were defined in paragraph 3.4.1 Since a and 
b are constant for x < Xi and for a; > x(cr^), it suffice to prove the estimates on 
{xi,x{a~). On this interval, a and 6 are defined through formulas of the type 

a{x{s))^a{s), b{x{s))^ I3{s). 

Differentiating these inequalities with respect to s, we obtain 

x'is)a'{x{s)) = a'is), 

and thus, using the formula ( A.l[ ), 

Iterating this process and using the fact that sin remains bounded away from zero 
in the interval under consideration, we infer that for all k G {0, • • • 4}, 

A; 
1=0 

and eventually, for all p G [1, cxd], 

k 
1=0 

Of course the same type of inequality holds for b as well. Therefore we now compute 
the derivatives of a and (3 with respect to s up to order 4. Since 

a{s) = -(1 - ^(s))V'?|ao(s) - mivis) - yi). 

we start with the derivatives of (3. The computations are lengthy but do not raise 
any difficulty. It can be easily checked that the most singular estimates correspond 
to the case when cos 9 vanishes algebraically near s ~ inf /i (assumption (H2i)), 
with the lowest possible exponent n. Indeed, if cos vanishes exponentially near 
s — inf Ji , then for all A; > 

||a,Vlk-=o(|inH''), 

while in case (i) 

||a^>IU~=o(z.-^). 

Similar formulas hold for S^V't'iao, ^s^n'/'t'ian- 
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Therefore in the rest of the proof, we only treat the case (H2i) with n = A 
(Remember that because of assumption (H4), cos, 9 vanishes at least like (s— inf /i)**) 
near inf /i). We explain with some detail the estimates of /3 and /?', and we leave 
the rest of the derivatives to the reader. 

First, we recall that by definition of a^ (see (3.37)), 



for all s e (cr , si). Therefore, using the estimates of Lemma 3.1.2 we have 



s) < 



so that 

V^"\moo^0{\\nv\-^)^0{1). 

Differentiating /? with respect to s, we obtain, for s G (c^, si), 
13' [s] = (!-¥') -d' cos 6'a„V°|aa - sin6ia^9„?/'?|an 
-0'sin09,^Oaf,+cos0a,V?|an 
+V' sin fiia,, Vt%o - 2 cos 6*9, V°|ao + ^' sin 6'V't°| an 
-93" cos 61 V'?! an > 



so that 

|^'(s)| < C(l-^) 

Eventually, we retrieve 



5y\\u5y\ 



5y 



5l\\n5y\ 



C\ip"\ |cos6i| 



mW=0{v-^''\\nv\-^). 
The same kind of estimate yields, for fc = 2, 3, 4, 

Now, since 

*= / jk-l 



1=0 



aW =^cH^(')ar^|ao- /?'')- 



(?/(g)-2yi) 



we infer, for fc > 1, 

lla^lUi < C||^W|Ui 



/=o*^-i 






-4+ 19 |lni/| 



ds 



fc-; 



L=°(rT-,Si 



92 



4. PROOF OF CONVERGENCE 



It can be checked that for / e {1, • • • , 4}, 

\\di^\dn\\LHs»,s+) < C\lniy\-^iy^, 
dKy{s)-yi) 



ds^ 



Q-l 



Gathering all the terms, wc obtain, for fc > 1, 

||aW|Ui=0(i/^). 
The estimates of a" and P" in L^ go along the same lines. 
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Notations 

Ai, Bi subdomains of 51 where the solution to the Sverdrup equation is continuous, 

v-mM 

r_E, Tw East and West parts of the boundary, p. [8] 

Fjv, Ts North and South parts of the boundary, p. [8] 

7i^i_i truncation of the boundary conditions to be lifted by the North/South bound- 
ary layers, p. [51] 



59 



li,i-i ^ li.i-i truncation of the North/South boundary layers, p 

5y — v^l^\ lni'p^(ln I Ini^l)^'' truncation parameter with respect to y, p 

5^' = XEiVi ± Sy) — xsiyi) local truncation parameter with respect to x, p 



41 



42 



9{s) angle between the horizontal vector Cx and the exterior normal to the boundary, 

p-Hl 

/+ set of indices corresponding to East corners {si € ^r^), p. [s] 

Ab, Xw inverse sizes of the East and West boundary layers, p. [2l] 

Aat == A5 = As = i^^^' * inverse size of the horizontal boundary layers, p. 23 

A4 majorizing function, p. [42] 

Pi partition of unity, p. [8[ 

s curvilinear abcissa, p. |8] 

Si curvilinear abcissa of the singular points {si G OTe U dVw), P- p^ 

Sj extremal points of the domain of validity of East and West boundary layers (if 



relevant), p. 23 

tf curvihnear abcissa of the East end of S, p. 33 
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CTj extremal points where energy is injected in North and South boundary layers 



(if relevant), p. 48 49 



S surfaces of discontinuity (ordonnates of critical points), p. 28 

T forcing term (coming from Ekman pumping), p. fl] 

iff truncation dealing with the transition between East/West and North/South 
boundary layer terms, p. Hs) |49] 



Xo localization near the boundary, p. [8] 
X± localization on f7*, p. 29 



Xv truncation of r near the singular point s^, p. 41 



a;^ = x^{y) graphs of the East boundaries of fl^, p. 
xY abcissa of the West end of S, p. 61 

(Xi 



29 



,yi) coordinates of the singular points {{xi,yi) = {x{si),y{si))), p. M 



ip^ solution to the transport equation, p. 40 



ipf solution to the transport equation with truncated source term, p. 



41 



ipM,s North and South boundary layer terms, p. 24 
'4'E,w East and West boundary layer terms, p. 22 



"0^'^* interior singular terms lifting the discontinuities of ip^, p. 
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ip'^ interior singular layer terms, p. 33 

V'int = "0? + "0^'^' regularization of the solution to the Sverdrup equation, p. 32 






ijj = ipf + TpE + V'b"' sum of all terms whose trace is lifted by the North/South 



boundary layers, p. 48 



ipn = V"''^' + V-'^X ( ^ s^^ I + V'w'"'' sum of all singular layer terms whose trace is lifted 



on the West boundary, p. 66 



Y parametrization of the interior singular layer, p. 33 
il = r^i \ ul^j^j domain with islands Hi, p 



11 



il subdomains of D, where the transport equation has a continuous solution, p. 28 



distance to the boundary, p. |8] 



Z — \{s)z scaled distance to the boundary, p. 20 
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